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1. Introduction 



Quantum physics is largely governed by the Schroedinger equation. Yet, exact solutions of 
the equation are relatively few. Besides lattice and other numerical calculations, we rely 
mostly on perturbative expansions. Such expansion quite often leads to a divergent series 
with zero radius of convergence, as in quantum electrodynamics, quantum chromodynamics 
and problems involving tunnelling and instantons. In a series of previous papers[l-4] we have 
presented a new approach to solve the low lying states of the Schroedinger equation. In 
the special case of one dimensional problems, this new approach leads to explicit convergent 
iterative solutions, in contrast to perturbative series expansions. These convergent solutions 
include the long standing difficult problem [5-14] of a quartic potential with symmetric minima. 

In this paper, we discuss some additional results bearing on the new method. In the one- 
dimensional case, we show that by changing the boundary condition to be applied at each 
iteration, we can obtain a convergent alternating sequence for the groundstate energy and 
wave function instead of the monotonic sequence found before[4]. This result will be spelled 
out later in this section and proved in Section 3. We also find that the asymmetric quartic 
double-well potential can be treated by an extension of the procedure used previously for the 
symmetric case. This extension is treated in Section 4. 

In addition, we have begun the exploration of higher dimensional problems along the same 
line. Although the same kind of iterative procedure can be set up, the linear inhomogeneous 
equation to be solved at each step cannot now be reduced to simple quadratures, as was 
done for one dimension. However, it is of interest that this equation is identical in form to 
an electrostatic analog problem with a given position dependent dielectric constant media; 
at each nth iteration, there is an external electrostatic charge distribution determined by the 
(n — l)th iterated solution, as we shall discuss in this section. 

Consider the Schroedinger equation 

Htp = Etl) (1.1) 

where H is the Hamiltonian operator, ip the wave function and E its energy. For different 
physics problems, H assumes different forms. For example, for a system of n non-relativistic 
particles in three dimensions, H may be written as 

H = J2C ijPiPj + V(x) (1.2) 

where x stands for xi, x 2 , • • • , x 3n the coordinate components of these n particles, V(x) is 
the potential function, dj are constants and p±, p 2 , Pz n are the momentum operators 
satisfying the commutation relation 

[pi, Xj] = -iSij. (1.3) 

(Throughout the paper, we set Planck's constant h = 1.) For a relativistic field theory, the 
Hamiltonian usually takes on a different form. Let $(r) be a scalar boson field at a three- 
dimensional position vector r, and n(r) be the corresponding conjugate momentum operator. 
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In this case we may write 

H = J d 3 r[U 2 {r) + V(Q(t))] (1.4) 
with II(r) and $(r') satisfying the commutation relation 

[n(r), $(r')] = -i5 3 (r-r'). (1.5) 

In both cases, the dependence of H on the momentum operators p t and n(r) are quadratic. 
Consequently, they can be brought into an identical standard form. In the above case of a 
system of non-relativistic particles, through a linear transformation 

{xi} - {ft}, (1.6) 
the Hamiltonian (1.2) can be written in the standard form 

H = ~V 2 + V( qi , q 2 , ftv) (1.7) 

with 

N Q2 

v 2 = E#2- (1-8) 

Likewise for the relativistic boson field Hamiltonian (1.4), we can use the Fourier-components 
of $(r) and n(r) as the set {x^ and {pi}. Through a similar transformation (1.6), the field 
Hamiltonian (1.4) can also be brought into the standard form (1.7) and (1.8), but with the 
number of variables N = oo. All our subsequent discussions will start from the Schroedinger 
equation in this standard form (1.7) - (1.8). Furthermore, in this paper, we shall limit our 
discussions only to the groundstate. 
In order to solve 

(_I V 2 + V(q)Mq) = E^q) (1.9) 
where q stands for the set {ft}, we proceed as follows: 

1. Construct a good trial function 0(q). A rather efficient way to find such trial functions is 
given in the next section. 

2. By differentiating 0, we define 

tf(q)-^M(q) _1 (ivV(q)), (1-10) 

in which the constant E may be determined by, e.g., setting the minimum value of U(q) to 
be zero. Thus, 0(q) satisfies a different Schroedinger equation 

(-^V 2 + f/(q))0(q)=£ o 0(q). (1.11) 
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Define w(q) and £ by 

U(q) = V(q)+w(q) (1.12) 

and 

E = E + £. (1.13) 
The original Schroedinger equation (1.9) can then be written as 

{-\V 2 + U(q) - £oMq) = Mq) - 5)^(q). (1.14) 

Multiplying this equation on the left by 0(q) and (1.11) by ip(q), we find their difference to 
be 

-«V • (0V^ - ^V0) = (w - £)V0- (1.15) 
The integration of its lefthand side over all space is zero, which yields 

£ --&f' (116) 

3. The above equation (1.14) will be solved iteratively by considering the sequences 

iki, fa, •", ikn, ■■■ and £ u £ 2 , ••-,£*, ••• (1-17) 

that satisfy 

(-^V 2 + C/(q) - So)Vn(q) = Mq) - W„-i(q) (1-18) 

with 

^o(q)=0(q). (1-19) 
As in (1.15) and (1.16), we multiply (1.11) by ip n and (1.18) by 0; their difference gives 

_I V • (0V^n - VW0) = (W - £n)^n-10 (1-20) 



and therefore 



(1.21) 



As we shall show, for many interesting problems 

lim £ n = £ and lim ib n = ib, (1-22) 



in contrast to the perturbative series expansion using w(q) as the perturbation. The key 
difference lies in the above expression (1.21) of £ n , which is a ratio, with both its numerator 
and denominator depending on the (n — l) th iterative solution ipn-i- 

4. There exists a simple electrostatic analog problem for the iterative equation (1.18). Assum- 
ing that ■0 n _i(q) has already been solved, we can determine E n through (1.21). The righthand 
side of (1.20), defined by 

<7„(q) = (w(q) - £ n )^„-i(q)0(q), (1.23) 
is then a known function. Introduce 

/n(q)=^n(q)Mq)- (1-24) 
In terms of / n (q), the n th order iterative equation (1.20) becomes 

-^V • (0 2 V/ n ) = a n . (1.25) 

Consider a dielectric medium with a dielectric constant dependent on q, given by 

/<q)=0 2 (q). (1.26) 

Interpret <r n (q) as the external electrostatic charge distribution, |/„ the electrostatic potential, 
— |V/ n the electrostatic field and 

D n = ~KVf n (1.27) 

the corresponding displacement vector field. Thus (1.25) becomes 

V-D n = a n , (1.28) 

the Maxwell equation for this electrostatic analog problem. 

At infinity, <f>(oo) = 0. In accordance with (1.26) - (1-27), we also have D n (oo) = 0. 
Hence the integration of (1.28) leads to the total external electrostatic charge to be also zero; 
i.e., 

J <r n (q) d N q = (1.29) 

which is the same result given by (1.21) for the determination of £ n . Because the dielectric 
constant /c(q) in this analog problem is zero at q = oo, the dielectric media becomes a perfect 
dia-electric at oo. Thus, the equation of zero total charge, given by (1.29), may serve as a 
much simplified model of charge confinement, analogous to color confinement in quantum 
chromodynamics. 

We note that (1.25) can also be derived from a minimal principle by defining 

J(/ B (q)) = J {^(V/„) 2 + a n f n }d N q. (1.30) 
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Because of (1.29), the functional 7(/ n (q)) is invariant under 

/n(q) -> /n(q) + constant. (1.31) 

Since the quadratic part of i"(/ n (q)) is the integral of the positive definite |k(V/„) 2 , the 
curvature of 7(/„(q)) in the functional space / n (q) is always positive. Hence, 7(/ n (q)) has 
a minimum, and that minimum determines a unique electrostatic field — |V/ n , as we shall 
see. To establish the uniqueness, let us assume two different V/ n , both satisfy (1.25), with 
the same k = 2 and the same a n ; their difference would then satisfy (1.25) with a zero 
external charge distribution. For a n = 0, the minimum of 7(/„(q)) is clearly zero with the 
corresponding V/ n = 0. To derive /„(q) from V/„, there remains an additive constant at 
each iteration. As we shall show, this arbitrariness allows us the freedom to derive different 
types of convergent series. 

To illustrate this freedom, let us consider a one-dimensional problem in which we may 
replace the variables {qi} by a single x. Furthermore, for this discussion, let us assume the 
potential V(x) to be an even function, with 

V(x) = V(-x) (1.32) 

(a condition that will be relaxed in our later analysis). The evenness of V(x) requires tp(x) = 
ip(—x) and therefore also 4>(x) = <fi(—x). Thus, we need only to consider the half-space 

x > 0. (1.33) 

Equations (1.24), (1.27) and (1.28) can be written now as 

Mx)=*l) n (x)/<i>(x), (1.34) 

Dn = -\K{x)f n {x) (1.35) 

and 

D' n {x)=a n {x) (1.36) 

where 

k{x) = (j) 2 (x) (1.37) 

and 

a n (x) = (w(x) - £ n )(i){x)i) n _i{x) (1.38) 
same as before. Throughout the paper, ' denote ^. 
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From (1.36) and D n (oo) = 0, we have 

D n (x) = - / a n (z)dz (1.39) 

Jx 

and, since a n (x) is even in x, we have from (1.29), 

roo 

/ a n (z)dz = 0. (1.40) 
Jo 

It follows then from (1.34) and (1.38)-(1.40), 

Dn(x) = -J^(w(z)-S n ) ( f> 2 (z)f n ^(z) dz 

and (1.41) 

D n{x) = Jo(w(z) - E n )(\> 2 {z)j n _ 1 {z) dz, 

which lead to, through (1.35), 

f n (x) = /„(oo) -2 r<p-\y) !°° 4> 2 {z){w{z)-£ n )f n ^{z)dz (1.42) 

Jx Jy 

and 

Ux) = /„(0) -2 r0- 2 (y) f V <p 2 {z){w{z)-£ n )f n ^{z)dz. (1.43) 

JO JO 

Consider first the case that in (1.38) is positive and satisfies 

w'(x) < for x>0. (1.44) 

The hierarchy theorem that will be proved in Section 3 states that if w(x) satisfies (1.44) then 
the iterative solution of (1.42) with the boundary condition 

/ n (oo) = 1 for all n (1.45) 

gives a convergent monotonic sequence £ 1; £ 2 , £3, ■ ■ •, where for all n, 

£n > £n-l, (1-46) 

and 

£ n — > £ as n^oo; (1-47) 

likewise, the sequence fo(x) = 1, /2(x), • • • is also monotonic and convergent at any 

x > with 

fn{x) > fn-l{x) (1.48) 



and 



f n (x) -> f(x) as n -> oo. (1.49) 

Furthermore, the convergence of (1.47) and (1.49) can hold for arbitrarily large but finite 
u>(x). A result that is surprising, but pleasant. 

On the other hand, if instead of (1.45), we impose a different boundary condition, one 
given by 

/ n (0) = 1 for all n, (1.50) 

then instead of (1.47), we have for all odd n = 2m + 1 an ascending sequence 

£i <8 3 <8 5 < (1.51) 
however, for the even n = 2m series, we have a descending sequence 

8 2 >8 4 >8 % > (1.52) 
furthermore, between any even n = 2m and any odd n — 21 + 1, we have 

8 2 m > &21+1 ■ (1.53) 

Since according to (1.13), 8 n is the n th order iteration towards 

8 = Eq-E, (1.54) 

each odd member £21+1 in (1.51) gives an upper bound of E, whereas each even member 8 2m 
in (1.52) leads to a lower bound of E. Both sequences approach the correct 8 as n — > 00, 
one from above and the other from below. For the boundary condition (1.50), our proof of 
convergence requires a condition on the magnitude of w(x). Still this is quite a remarkable 
result. 

In Section 2, we discuss the details of how to construct a good trial function 0(q) for 
the N-dimensional problem. Section 3 gives the proof of the hierarchy theorem for the one- 
dimensional problem in which V(x) = V(— x) is an even function of x and the potential- 
difference function w(x) is assumed to satisfy (1.44); i.e., w'(x) < for x > 0. The extension 
to the asymmetric case V(x) 7^ V(—x) is discussed in Section 4. The hierarchy theorem is 
also applicable to Mathieu's equation, which has infinite number of maxima and minima. In 
the Appendix, we give a soluble example in one dimension. 

In dimensions greater than 1, at each iteration Eq.(1.21) gives a fine tuning of the energy, 
just like the one-dimensional problem. Hence, there are good reasons to expect our approach 
to yield convergent solutions in any higher dimension. In Section 5, we formulate an explicit 
conjecture to this effect. We describe an attempt to prove this conjecture by generalizing the 
steps used to prove the hierarchy theorem in one dimension. The attempt fails at present 
because the proof of one of the lemmas does not appear to generalize in higher dimension. 

The present paper represents the synthesis and generalization of results, some of which 
have appeared in our earlier publications[l-4]. The function D n introduced in this paper is 
identical to the function h n used in Ref.[4]. 
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2. Construction of Trial Functions 



2.1 A New Formulation of Perturbative Expansion 

In many problems of interest, perturbative expansion lead to asymptotic series, which is 
not the aim of this paper. Nevertheless, the first few terms of such an expansion could provide 
important insight to what a good trial function might be. For our purpose, a particularly 
convenient way is to follow the method developed in Refs.fl] and [2]. As we shall see, in this 
new method to each order of the perturbation, the wave function is always expressible in terms 
of a single line-integral in the N-dimensional coordinate space, which can be readily used for 
the construction of the trial wave function. 

We begin with the Hamiltonian H in its standard form (1.7). Assume V(q) to be positive 
definite, and choose its minimum to be at q = 0, with 

V(q) > V(0) = 0. (2.1) 

Introduce a scale factor g 2 by writing 

V(q) = g 2 v(q) (2-2) 

and correspondingly 

^(q) = e -f^). (2.3) 
Thus, the Schroedinger equation (1.9) becomes 

(_Iy 2 + g 2 v(q))e~ 9S ^ = Ee^ s ^ (2.4) 
2 

where, as before, q denotes q 1: q 2) • • • , ?jv ar| d V the corresponding gradient operator. 
Hence S(q) satisfies 

- l -g 2 (VS) 2 + l -gV 2 S + g 2 v = E. (2.5) 
Considering the case of large g, we expand 

S(q) = 5 (q) + ^(q) + <T 2 S 2 (q) + • • • (2.6) 

and 

E = gE + E 1 + g~ 1 E 2 + ■■■. (2.7) 

Substituting (2.6) - (2.7) into (2.5) and equating the coefficients of g~ n on both sides, we 
find 

(VSo) 2 = 2v, 
VS -VS 1 = ^V 2 S -E , 

VS -VS 2 = ±[V 2 S 1 -(VS 1 ) 2 ]-E 1 , (2.8) 
VSo • VS 3 = \ [V 2 S 2 -2V5i- VS 2 ] - E 2 , 
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etc. In this way, the second order partial differential equation (2.5) is reduced to a series of 
first order partial differential equations (2.8). The first of this set of equations can be written 
as 

1 -[VS ( q )f-v( q ) = + . (2.9) 

As noted in Ref.[l], this is precisely the Hamilton-Jacobi equation of a single particle with 
unit mass moving in a potential "— i>(q)" in the N-dimensional q-space. Since q = is 
the maximum of the classical potential energy function — v (q), for any point q ^ there 
is always a classical trajectory with a total energy 0+, which begins from q = and ends 
at the other point q ^ 0, with S (q) given by the corresponding classical action integral. 
Furthermore, Sb(q) increases along the direction of the trajectory, which can be extended 
beyond the selected point q ^ 0, towards oo. At infinity, it is easy to see that S (q) = oo, 
and therefore the corresponding wave amplitude e~ 95 °( q ) is zero. To solve the second equation 
in (2.8), we note that, in accordance with (2.1) - (2.2) at q = 0, VS oc v*(Q) = 0. By 
requiring Si(q) to be analytic at q = 0, we determine 

£o = ^(V 2 S ) atq=0 . (2.10) 

It is convenient to consider the surface 

S (q) = constant; (2-11) 

its normal is along the corresponding classical trajectory passing through q. Characterize each 
classical trajectory by the S -value along the trajectory and a set of iV — 1 angular variables 

a = («i(q), a 2 (q), • • • , ajv-i(q)), (2.12) 
so that each a determines one classical trajectory with 

Va.rVSo = 0, (2.13) 

where 

j = l,2,---,iV-l. (2.14) 

(As an example, we note that as q — > 0, t>(q) — > and therefore S — > ^J2^iQi- 

i i 

Consider the ellipsoidal surface So = constant. For So sufficiently small, each classical trajec- 
tory is normal to this ellipsoidal surface. A convenient choice of a could be simply any N — 1 
orthogonal parametric coordinates on the surface.) Each a designates one classical trajectory, 
and vice versa. Every (So, a) is mapped into a unique set (q 1: q 2 , ■ ■ ■ , Qn) with S > by 
construction. In what follows, we regard the points in the q-space as specified by the coor- 
dinates (S , a). Depending on the problem, the mapping (q 1 , q 2 , gjv) — ► (So, a) may 
or may not be one-to-one. We note that, for q near 0, different trajectories emanating from 
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q = have to go along different directions, and therefore must associate with different a. 
Later on, as So increases each different trajectory retains its initially different a-designation; 
consequently, using (So, a) as the primary coordinates, different trajectories never cross each 
other. The trouble-some complications of trajectory-crossing in q-space is automatically re- 
solved by using (So, a) as coordinates. Keeping a fixed, the set of first order partial differential 
equation can be further reduced to a set of first order ordinary differential equation, which are 
readily solvable, as we shall see. 
Write 

S 1 (q)=S 1 (S ,a), (2.15) 

the second line of (2.8) becomes 

(VS ) 2 (^r) a = -V 2 S -E , (2.16) 
and leads to, besides(2.10), also 

5 1 (q) = 5 1 (5 ,a) = / J^i^So-Eo], (2.17) 

where the integration is taken along the classical trajectory of constant a. Likewise, the third, 
fourth and other lines of (2.8) lead to 

E 1 = h L V 2 S l - (V^) 2 ] at q=0 , (2.18) 

S 2 ( q ) = S 2 (S ,a) = J-^ 2 { 1 -[vz Sl -(VS i y}-E 1 }, (2.19) 

E 2 = ^[V 2 S 2 - 2(VS X ) • (VS 2 )] at q=0 , (2.20) 

Ss( q ) = S 3 (S ,a) = Jj^{-[V 2 S 2 -2(VS 1 )-(VS 2 )]-E 2 }, (2.21) 

etc. These solutions give the convenient normalization convention at q = 0, 

S(0) = 

and 

e- s(0) = 1. (2.22) 
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Remarks 



(i) As an example, consider an N-dimensional harmonic oscillator with 

V(q) = ^(q! + q 2 2 + --- + q 2 N ). (2.23) 

From (2.2), one sees that the Hamilton-Jacobi equation (2.9) is for a particle moving in a 
potential given by 

-v(cL) = -\(q 2 + q 2 2 + --- + q 2 N ). (2.24) 

Thus, for any point q^O the classical trajectory of interest is simply a straight line connecting 
the origin and the specific point, with the action 

S^) = \{ql + ql + --- + q 2 N )- (2-25) 
The corresponding energy is, in accordance with (2.10), 

Eo = y- (2.26) 

By using (2.8), one can readily show that E x — E 2 — ■ ■ ■ — and Si = S 2 = ■ ■ ■ = 0. The 
result is the well known exact answer with the groundstate wave function for the Schroedinger 
equation (2.4) given by 

e -«*(i> = exp [- 9 -{q\ + q 2 2 + --- + q 2 N )} (2.27) 
and the corresponding energy 

N 

E = -g. (2.28) 

(ii) From this example, it is clear that the above expression (2.6) - (2.8) is not the well-known 
WKB method. The new formalism uses — f(q) as the potential for the Hamilton-Jacobi 
equation, and its "classical" trajectory carries a 0+ energy; consequently, unlike the WKB 
method, there is no turning point along the classical trajectory, and the formalism is applicable 
to arbitrary dimensions. 
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2.2 Trial Function for the Quantum Double- well Potential 



To illustrate how to construct a trial function, consider the quartic potential in one dimen- 
sion with degenerate minima: 

V{x) = \g 2 {x 2 -a 2 ) 2 . (2.29) 

An alternative form of the same problem can be obtained by setting q = y/2ga(a — x) so that 
the Hamiltonian becomes 

H = ~\ ^ + V{X) = 29<Xk) (2 ' 30) 

where 

This shows that the dimensionless (small) expansion parameter is related to \/y/&gd 6 ; as it 
turns out, the relevant parameter is its square. In the following, we shall take a = 1 so that 
the expansion parameter is 1/g; in the literature[5-14] one often finds the assumption 2ga = 1 
(placing the second minimum of the potential at q — 1/g) so that l/y/8ga? reduces to g and 
the anharmonic potential appears as (l/2)g 2 (l — gq) 2 . Then g appears with positive powers 
instead of negative, but the coefficients of the power series are the same as with our form of 
the potential, apart from the overall factor 2ga. 

For the above potential (2.29), the Schroedinger equation (2.4) is (with a = 1) 

{ ~\^ + ~ l?)i){x) = m%) (2 ' 32) 

where, as before, ip(x) = e~ gS ^ is the groundstate wave function and E its energy. Using 
the expansions (2.6) - (2.7) and following the steps (2.8), (2.10) and (2.15) - (2.21), we find 
the well-known perturbative series 

S,{x)= l -{x-l)\x + 2l S 1 (x) = \n X ±±, S 2 (x) = 1 - ^^~ 2 , • • <2.33) 

and 

Eq = 1, E 1 = ~ ^2 = ~, •••• (2.34) 

Both expansions S = S + g~ 1 Si + g~ 2 S 2 -\ and E = gE + E 1 + g~ l E 2 -\ are divergent, 

furthermore, at x — — 1 and for n > 1, each S n (x) is infinite. The reflection x — > — x gives 
a corresponding asymptotic expansion S n (x) — > S n (—x), in which each S n (—x) is regular at 
x = —1, but singular at x = +1. 
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We note that for g large, the first few terms of the perturbative series (with (2.33) for 
x positive and the corresponding expansion S n (x) — > S n (—x) for x negative) give a fairly 
good description of the true wave function ip(x) whenever ip(x) is large (i.e. for x near ±1). 
However, for x near zero, when ij)(x) is exponentially small, the perturbative series becomes 
totally unreliable. This suggests the use of first few terms of the perturbative series for regions 
whenever ip(x) is expected to be large. In regions where tp(x) is exponentially small, simple 
interpolations by hand may already be adequate for a trial function, as we shall see. Since the 
quartic potential (2.29) is even in x, so is the groundstate wave function; likewise, we require 
the trial function <f>(x) also to satisfy 4>(x) = (/>(—x). At x = 0, we require 

(^t) x=0 = 0'(O) = 0. (2.35) 

To construct <f>(x), we start with the first two functions S (x) and Si(x) in (2.33). Introduce, 

for x > 0, 

<f) + (x) = e-ySoW-s^x) = (_^_) e -<7So(*) (2.36) 

and 

In order to satisfy (2.35), we define 

f (f)+(x) + a=±<i>Jx), for 0<x<l 
= = for x ~ >1 (2.38) 

Thus, by construct 0'(O) = 0, <p(x) is continuous everywhere, for x from — oo to oo, and so 
is its derivative. 

By differentiating <p + (x) and <p(x), we see that they satisfy 

(T + V + u)4> + = g(f> + (2.39) 

and 

(T + V + w)<f) = g<j), (2.40) 

where 
and 

w (x) = w(-x) (2.42) 



14 



with, for x > 

w(x) = u(x) + g(x) (2.43) 

where 

r (g _ 1)e29 s oW -f 9 

= ^ViHto-De^oM-f.' " (2.44) 

[ for x > 1. 

Note that for > 1, ^(x) is positive, and has a discontinuity at x = 1. Furthermore, for re 
positive both m(x) and ^(rr) are decreasing functions of x. Therefore, w(x) also satisfies for 

x > 0, 

iu'(a;) < 0, (2.45) 

a property that is very useful in our proof of convergence which will be discussed in the next 
section. 
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3. Hierarchy Theorem and Its Generalization 



In this section, we restrict our discussions to a one-dimensional problem, in which the potential 
V(x) is an even function of x, as in the example given in the previous Section 2.2. The 
Schroedinger equation (1.9) becomes 

-\^"{x) + V(x)if>(x) = E^{x) (3.1) 

with ip(x) as its groundstate wave function, E the groundstate energy and ' denoting ^, as 
before. For the one-dimensional problem, the trial function <f)(x) satisfies 

-±<j>»(x) + U(x)<i>(x) = Eo<i>(x), (3.2) 
as in (1.11); therefore (3.1) can be written as 

+ (U(x) - EoMx) = (w(x) - EMx), (3.3) 



in which 



and 



U(x)=V(x)+w(x) (3.4) 



E = E + S, (3.5) 
as before. Throughout this section, we assume 

V(x) = V(-x), U(x) = U(-x), 

ip(x) = il)(—x) and (j>(x) — (/>(— x); (3.6) 
hence, we need only to consider 

x > 0. (3.7) 

Furthermore, as in the example of the symmetric quartic double-well potential given in Sec- 
tion 2.2, we assume w(x) to satisfy 

w'(x) < for x > (3.8) 

and 

w(oo) = 0. (3.9) 
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Therefore, w(x) is positive for x positive. Otherwise, the shape of w(x) can be arbitrary. The 
Schroedinger equation (3.1) will be solved through the iterative steps (1.34) - (1.43), using 
the sequences 

£i, £2, •••,£„, ••• (3.10) 

for the energy difference £ = E — E, and the sequence 

fl(x), f 2 (x), f n (x), ••• (3.11) 

for the ratio f(x) = ip(x)/(/)(x) with, for n = 0, 

f (x) = 1. (3.12) 

In this section, we differentiate two sets of sequences, labelled A and B, satisfying different 
boundary conditions: 

/ re (oo) = 1 for all n, in Case (A) 

or 

/„(0) = 1 for all n. in Case (B) 

Thus, in accordance with (1.42)-(1.43), we have in Case (A) 

f n (x) = 1 - 2 / r 2 (y)dy / 4>\z){w{z) - S n )f n ^(z)dz, (3.13A) 

Jx Jy 

whereas in Case (B) 

f n (x) = 1- 2 r<p-\y)dy f V <t> 2 {z){w{z)-£ n )f n _ 1 {z)dz. (3.135) 

JO JO 

In both cases, S n is determined by the corresponding f n -i(x) through (1.38) and (1.40); i.e., 

€ n = [w /„_i]/[/„_i] (3.14) 
in which [F] of any function F(x) is defined to be 

[F] = / <p 2 (x)F(x)dx. (3.15) 

JO 

Eqs. (1.35) and (1.37) give 

f n {x) = -2<p- 2 D n (x). (3.16) 

Likewise, (1.38) - (1.39) lead to 

[•CO 

D n (x) = - / <P 2 (z)(w(z) - £ n )j n _ x {z)dz (3.17) 

17 



which, on account of (1.40), is identical to 

D n (x) = f ( f> 2 (z)(w(z)-£ n )f n _ 1 (z)dz. (3.18) 
Jo 

These two expressions of D n (x) are valid for both cases (A) and (B). Let x n be defined by 

w(x) — £ n = at x = x n . (3.19) 

Since w'(x) < 0, (3.19) has one and only one solution, with w(x) — £ n negative for x > x n 
and positive for x < x n . Thus, if 

fn-i(x) > (3.20) 
for all x > 0, we have from (3.17) - (3.18) 

D n (x) > (3.21) 

and therefore, on account of (3.16), 

f n {x) < 0. (3.22) 

In terms of the electrostatic analog introduced in Section 1, through (1.26) - (1.29), one can 
form a simple physical picture of these expressions. Represent D n (x) by the standard flux 
of lines of force. Because the dielectric constant k(x) = 4> 2 (x) is zero at x = oo, so is the 
displacement field. Hence, D n (po) = 0; therefore each line of force must terminate at a finite 
point. Since the electric charge density is a n (x) = 4> 2 (x)(w(x) —£ n )f n -i(x), the total electric 
charge to the right of x is 

Q n {x) = / a n (z)dz. 

J x 

It must also be the negative of the flux D n (x) passing through the same point x: i.e., 

POO 

D n (x) = -Q n (x) = - cr n (z)dz, 

Jx 

which gives (3.17). In the whole range from x = to oo, the total electric charge / °° a n (z)dz 
is zero; therefore, we have 

Q n (0) = and D n (0) = 0. 
Furthermore, at any point x > 0, the total charge from the origin to the point x is 

r-x 

\ a n (z)dz, 
Jo 

which must also be the negative of the above Q n (x), and therefore the same as D n (x); that 
leads to (3.18). From (3.19) and w'(x) < 0, one sees that the charge distribution a n (x) is 



18 



negative for x > x n , at x = x n and positive for < x < x n . Correspondingly, moving from 
x = oo towards the left, the displacement field increases from D n (oo) = to D n (x) > 0, 
reaching its maximum at x — x n , then as x further decreases, so does D n (x), and finally 
reaches D n (0) = at x = 0. 

In Case (A), because of f n (oo) = 1, (3.22) leads to 

/n(0) > f n (x) > / n (00) = I. (3.23A) 

Since for n — 0, /o(^) = 1> (3.20) - (3.23A) are valid for n = 1; by induction these expressions 
also hold for all n; in Case (A), their validity imposes no restriction on the magnitude of w(x). 
In Case (B) we assume w(x) to be not too large, so that (3.13B) is consistent with 

f n (x) > for i>0 

and therefore 



/n(0) = 1 > f n (x) > / n (00) > 0. (3.235) 

As we shall see, these two boundary conditions (A) and (B) produce sequences that have 
very different behavior. Yet, they also share a number of common properties. 

Hierarchy Theorem (A) With the boundary condition /„(oo) = 1, we have for all n 

£n+l > £n (3.24) 

and 

d ( Wx)\ < Q at any x > Q 
ax \ f n {x) J 

Thus, the sequences {£„} and {f n (x)} are all monotonic, with 

S 1 <S 2 <S 3 <--- (3.26) 

and 



l<h(x)<f 2 (x)<f 3 (x)<--- (3.27) 

at all finite and positive x. 

(B) With the boundary condition f n (0) = 1, we have for all odd n = 2m + 1 an ascending 
sequence 

S l < £ 3 < £ 5 < ■ ■ ■ , (3.28) 
but for all even n = 2m, a descending sequence 

S 2 > S A > S Q > ■■■■ (3.29) 
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furthermore, between any even n = 2m and any odd n — 21 + 1 



S 2m > Sn+i- (3.30) 
Likewise, at any x > 0, for any even n = 2m 

d ( f2m+i(xy 



whereas for any odd n = 21 + 1 



dx \ /2m (x) 
d / f2l+2( X ) 



< 0, (3.31) 



dx \f 2 i+i(x) 



> 0. (3.32) 



Remarks. 

1. The validity of Eqs. (3.24) and (3.25) for the boundary condition / n (oo) = 1 was estab- 
lished in Ref.[4]. The validity of Eqs. (3.28)-(3.32) for the boundary condition / n (0) = 1 is 
the new result of this paper, which we shall establish. 

2. As we shall also show, the lowest eigenvalue E of the Hamiltonian T + V is the limit of 
the sequence {E n } with 

E n = E — £ n . (3.33) 
Thus, the boundary condition / n (oo) = 1 yields a sequence, in accordance with (3.26), 

E 1 > E 2 > E 3 > ■ ■> E, (3.34) 

with each member E n an upper bound of E, similar to the usual variational iterative sequence. 

3. On the other hand, with the boundary condition / n (0) = 1, while the sequence of its odd 
members n = 21 + 1 yields a similar one, like (3.34), with 

Ei > E 3 > E 5 > ■ ■ ■ > E, (3.35) 

its even members n = 2m satisfy 

E 2 <E 4 <E 6 <--- <E. (3.36) 

It is unusual to have an iterative sequence of lower bounds of the eigenvalue E. Together, 
these sequences may be quite efficient to pinpoint the limiting E. 

The proof of the above generalized hierarchy theorem depends on several lemmas that are 
applicable to both boundary conditions: (A) /„(oo) = 1 and (B) / n (0) = 1; these lemmas 
will be established first, and then followed by the proof of the theorem. 
Lemma 1 For any pair f m (x) and fi(x) 

i) if at all x > 0, ^ (j^y) < then Sm+1 > ( 3 - 37 ) 
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and 



d ( f (x)\ 

ii) if at all x > 0, — > then E ™+i < £i+i- (3.38) 



Proof 



According to (3.14) 



Also by definition (3.15), 



Their difference gives 



£m+l[fm] = [W f m ]. (3.39) 



£l+l[fm] = [£l+l fm}- (3.40) 



(£ m+1 - £ l+1 )[f m ] = [(w - £ l+1 ) f m }. (3.41) 

From (3.14), 

= [{w- £ l+1 ) /,]. (3.42) 

Let xi + i be defined by (3.19). Multiplying (3.41) by fi(xi +1 ) and (3.42) by f m (xi +1 ) and 
taking their difference, we have 

fi(xi +1 )(£ m+1 - £i+i)[f m ] = [(w - £ l+1 ) (f m (x)fi(xi +1 ) - fi(x)f m (x l+1 ))], (3.43) 

in which the unsubscripted x acts as a dummy variable; thus [f m (x)] means [f m ] and [f m (xi +1 )] 
means f m (x l+1 ) • [1], etc. 

(i) If (f m (x)/fi(x))' < 0, then for x < x l+1 

fm{x) fm{xi+l) . , 

fi(x) fi(x l+1 ) ■ {6 ' } 

In addition, since w'(x) < and w(xi + ±) = £i+\, we also have for x < xi + i 

w(x) > £ l+1 . (3.45) 

Thus, the function inside the square bracket on the right hand side of (3.43) is positive 
for x < xi + \. Also, the inequalities (3.44) and (3.45) both reverse their signs for x > xi + \. 
Consequently, the right hand side of (3.43) is positive definite, and so is its left side. Therefore, 
on account of (3.23A)-(3.23B), (3.37) holds. 

(ii) If (fm(x)/fi(x))' > 0, we see that for x < xi+i, (3.44) reverses its sign but not (3.45). 
A similar reversal of sign happens for x > xi + \. Thus, the right hand side of (3.43) is now 
negative definite and therefore £ m +i < £i+i- Lemma 1 is proved. 
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The following lemma was already proved in Ref.[4]. For the convenience of the readers, 
we also include it in this paper. Let 

77 = 77(0 (3.46) 
be a single valued differentiable function of £ in the range between a and b with 

< a < £ < b (3.47) 

and with 

77(a) > 0. (3.48) 

Lemma 2 . 

(i) The ratio 77/^ is a decreasing function of £ for a < £ < b if 

at £ = a (3.49) 

and 

— '- < for a < £ < b. (3.50) 

(ii) The ratio 77/^ is an increasing function of £ for a < £ < b if 

|>! at { = . (3.51) 



and 

cPr) 

Proof Define 



to be the Legendre transform We have 



2 > for a < £ < b. (3.52) 



L = ^~V (3.53) 



and 



d 4 = &. (3.54) 



(3.55) 



Since (3.49) says that L(a) < and (3.50) says that J < for a < £ < 6, these two 
conditions imply L(£) < for a < £ < b, which proves (i) in view of (3.55). The proof of (ii) 
is the same, but with inequalities reversed. 
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Lemma 3 For any pair f m (x) and fi(x) 
(i) if over all x > 0, 

and (ii) if over all x > 0, 

Proof From (3.17)-(3.18), we have 

^m+lW = M&) -^m+l) 2 (^) /raW 

and 

D' l+1 (x) = (w(x) - 2 (x) /,(*). 

Define 

f = A+i(aO and 77 = D m+ i(x). 

In any local region of x where D' l+1 (x) 7^ 0, we can regard rj = r}(£) through rj(x) 
Hence, we have 

dv _ D' m+1 (x) =r(x) f m (x) 



where 



and 



di D[ +1 {x) v ' fi(x) 

w(x) - £ m +i 
w(x) - £ l+1 

±(M 1 (D'm + l\_ 1 (fm\ 

^U/ " D[ +1 \D[ +1 ) D[ +1 \ fj 

1 I 1 fm , ( fm \ \ 

■ /y _l_ /y ' 11 



D' l+1 \ fi \fi 



where 



£to+i — £1+1 _j, 
-1+1) 



// \ ■-'7711-1 ■-'(1-1 lj \ 

r (X) — ~, — / — \ ~ TTrW (x). 
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(i) If (fm/ fi)' < 0, from Lemma 1, we have 

S m+ i > £ i+1 . (3.65) 
From w'(x) < and the definition of x m+ \ and xi + ±, given by (3.19), we have 

x m+1 < x i+1 , (3.66) 

w{x m+ i) = £ m+1 and w(x i+1 ) = £ l+1 . (3.67) 

We note that from (3.17) - (3.18) D m+ i(x) and D i+ i(x) are both positive continuous 
functions of x, varying from at x — 0, 

D m+1 (0) = A+i(0) = (3.68) 

to at x = oo 

An + i(oo) = A+i(oo) = (3.69) 

with their maxima at x m+1 for D m+1 (x) and x l+1 for since in accordance with 

(3.58)-(3.59) and (3.67), 

D' m+1 {x m+l ) = and D' l+1 (x l+1 ) = 0. (3.70) 

From (3.64)-(3.65), we see that r'(x) is always < 0. Furthermore, from (3.62), we also find 
that the function r(x) has a discontinuity at x = xi + ±. At x = 0, r(0) satisfies 

< .(0) = W ^ - £ ; +1 < 1. (3.71) 
w(0) - 

As x increases from 0, r(x) decreases from r(0), through 

r{x m+1 ) = 0, (3.72) 

to — oo at x = xi + i— ; r(x) then switches to +oo at x = xi +1 + , and continues to decrease 
as x increases from xi + i+. At x = oo, r(x) becomes 

r(oo) = §i±i > 1. (3.73) 

It is convenient to divide the positive rr-axis into three regions: 

(I) < x < x m+1 , 

(II) x m+ i < x < xi+i (3.74) 

(III) x l+l < x. 
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region 




D' l+1 (x) 


w(x) - £ m+1 


w(x) - £ l+1 


r(x) 


r'(x) 


I 


> 


> 


> 


> 


> 


< 


II 


< 


> 


< 


> 


< 


< 


III 


< 


< 


< 


< 


> 


< 



Table 1. The signs of D' m+l (x), D' l+1 ( x) , w{x^ — £ m +\, w{x^ — £i+i, 
r(x) and r'(x) in the three regions defined by (3.74), when £ m+ i > £i + \. 



Table 1 summarizes the signs of D' m+1 , D' l+1 , r and r' in these regions. Assuming 
(fm/fi)' < we shall show separately the validity of (3.56), (An+i/A+i)' < 0, in each 
of these three regions. 

Since 

£j+i < w(x) < £ m+ i in II, (3.75) 

D m+ i(x) is decreasing and D t+1 (x) is increasing; it is clear that (3.56) holds in II. 

In each of regions (I) and (III), we have r(x) > from (3.62) and r'(x) < from 
(3.64). Since (f m / fi)' is always negative by the assumption in (3.56), both terms inside the 
big parenthesis of (3.63) are negative; hence the same (3.63) states that d 2 rj/d^ 2 has the 
opposite sign from D' l+1 . From the sign of D' l+1 listed in Table 1, we see that 



2 < in (I) (3.76) 



2 > in (III). (3.77) 



d 2 r] 

and 

d 2 r) 

Within each region, rj = D m+1 (x) and £ = Di+i(x) are both monotonic in x; therefore, rj is 
a single-valued function of £ and we can apply Lemma 2. In (I), at x = 0, both D m+1 (0) and 
D l+1 (0) are according to (3.18), but their ratio is given by 

An+i(Q) = DL+M ,o 7H) 

A + i(o) AVi(o) ' 1 ' ) 

Therefore, 

=(?) • (3.79) 

Furthermore, from (3.76), ^ < in (I), it follows from Lemma 2, case(i), the ratio r]/£_ is a 
decreasing function of £. Since £' = is > in (I), according to (3.59), we have 

<0 in (I).. (3.80) 
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In (III), at x = oo, both D rn+ i(oo) and Z^ +1 (oo) are according to (3.69). Their ratio is 

An+l(oo) _ D' m+1 (oo) 



which gives at x = oo, 



A+i(oo) 



dr]^ 



AVi(oo) ' 



?7 



(3.81) 



As x decreases from x = 00 to x = xi+i, from (3.77) we have ^ > in (III). It follows 



from Lemma 2, case (ii), 77/^ is an increasing function of £. Since £' = h' l+1 is < 0, because 
x > xi+i, we have 



d_ ( Drn+l 

dx \ A+i , 



< in (III).. 



(3.82) 



Thus, we prove case(i) of Lemma 3. Case(ii) of Lemma 3 follows from case (i) through the 

interchange of the subscripts m and I. Lemma 3 is then established. 

Lemma 4 Take any pair f m (x) and f t (x) 

(A) For the boundary condition f n (oo) = 1, if at all x > 0, 



d (fm(x)\ d ( f m+ l{x)' 

1 < then at all x > 0, 



dx \ fi(x) 
therefore, if at all x > 0, 



dx \ fi+i(x) 



<0; 



d ( fm(x) \ Q thgn ata \\ x> Q d ( fm+l( X Y 

dx V fi(x) J dx \ fi+i(x) 

(B) For the boundary condition /„(0) = 1, if at all x > 0, 



> 0. 



d_ ( fm(x) 

dx V fi(x) 
therefore, if at all x > 0, 

d_ ( fmjx) 

dx \ fi(x) 
Proof Define 

From (1.35) we see that 



< then at all x > 0, — 



d ( f m +i(x) 



> then at all x > 0, 



dx \ fi+i(x) 

d_ ( fm+l(x) 

dx V fi+i(x) 



>0; 



< 0. 



i = fi+i(x) and r) = f m+1 (x). 



drj _ f' m+1 {x) _ D m+1 (x) 



di fl +1 (x) D l+1 (x) 
26 



(3.83A) 



(3.8471) 



(3.83B) 



(3.845) 



(3.85) 



(3.86) 



and 

d ( df}\ 1 d f D m+ i(x) 



di\diJ fl +1 dx\D l+1 i 



x 



(3.87) 



(A) In this case. f n (oo) = 1 for all n. Thus, at x = oo, £ = /i +1 (oo) = 1, 17 = / m+1 (oo) = 1, 
and their ratio 

(I) = 1. (3.88) 

\ S / x=oo 

At the same point x = 00, in accordance with (3.17), Z? J+1 (oo) = D m+1 (oo) = 0, but their 
ratio is, on account of w(oo) = and (3.37) of Lemma 1, 



D m+1 (oo) _ D' m+1 (oo) _ w(oo) 

A+i(oo) D'i+i(oo) w(oo)-£ i+1 £ l+1 



> 1, (3.89) 



in which the last inequality follows from the same assumption, if (f m / //)' < 0, shared by 

(3.37) of Lemma 1 and the present (3.83A) that we wish to prove. Thus, from (3.86), at 

x = 00 

(I) =§^(1) = * (3 90 > 

As x decreases from 00 to 0, £ increases from // +1 (oo) = 1 to /;+i(0) > 1, in accordance with 
(3.22) and (3.23A). On account of (3.56) of Lemma 3, we have (D m+1 /D l+1 )' < 0, which 
when combined with (3.87) and /^(x) < leads to 

l( d 4)>0. (3.91) 
Thus, by using (3.51)-(3.52) of Lemma 2, we have rj/£ to be an increasing function of £; i.e., 

|(|) > a 



Because 



and f[ +1 < 0, we find 



=(fe)-=(!H 
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which establishes (3.83A). Through the interchange of the subscripts m and I, we also obtain 
(3.84A). 

Next, we turn to Case (B) with the boundary condition / n (0) = 1 for all n. Therefore at 

x = 0, 

**g = l. ,3,5) 

Furthermore from (3.16) and (3.18B), we also have f' n+1 (0) = // +1 (0) = and D m+1 (0) = 
D l+1 (0) = 0, with their ratio given by 



df m+ i\ _ /; t+ i(0) _ An+i(0) _ D' m+1 (0) 

(3.96) 



dfi + i) x=0 fl + M " A+i(0) " AVi(O) 

w(o) - £ m+i 



^(0) - s l+1 ' 

From (3.37) of Lemma 1, we see that if (fm/fi)' < 0, then £ m +i > £i+i and therefore 

< 1. (3.97) 

V d Ji+i ) x=0 

= 1. (3.98) 



x=0 



Thus, 

2=0 V£/ x=0 



^ <(i) . (3.99) 



Analogously to (3.53), define 

'; - (-.1. l •> ■ ,,, 

di fl +1 (x) 

therefore 



L(x) =i^L-fi = f l+ i(x)^±44 - f m+1 (x); (3.100) 



dL(x) £ t dL t,t d dr) 
dx ~ ~dj~ ~dj~dj 

= ^ (I) = fl+l i (^f) = fi+i i (^r) • (3 - 101) 

From (3.56) of Lemma 3, we know that if (f m /fi)' < then (D rn+ i/Di +1 )' < 0, which leads 
to 

^P- < 0. (3.102) 
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From (3.100), we have 

and therefore at x = 0, because of (3.99), 

L(0) < 0. (3.104) 
Combining (3.102) and (3.104), we derive 

L(x) < for i>0. (3.105) 
Multiplying (3.100) by f l+1 {x), we have 

f' l+1 (x)L{x) = f l+l (x)f m+1 (x) - f m+1 (x)fi +1 (x) 

Because f[ +1 (x) and L(x) are both negative, it follows then 

/ f m +i{x)\ 



>0, 



which gives (3.83B) for Case (B), with the boundary condition f n (0) = 1. Interchanging the 
subscripts m and I, (3.83B) becomes (3.84B), and Lemma 4 is established. 
We now turn to the proof of the theorem stated in (3.24)-(3.32). 

Proof of the Hierarchy Theorem 
When n = 0, we have 

f (x) = 1. (3.107) 

From (3.20)-(3.22), we find for n = 1 

f[{x) < 0, (3.108) 

and therefore 

(fi/fo)' < 0. (3.109) 

In Case (A), by using (3.83A) and by setting m = 1 and / = 0, we derive (/2//1)' < 0; 
through induction, it follows then (f n +i/ fn)' < for all n. From Lemma 1, we also find 
S n+ i > £ n for all n. Thus, (3.24)-(3.27) are established. 
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In Case (B), by using (3.109) and (3.83B), and setting m — 1 and / = 0, we find 
(iV/i)' > 0, which in turn leads to (/ 3 // 2 )' < 0, • • •, and (3.31)-(3.32). Inequalities (3.28)- 
(3.30) now follow from (3.37)-(3.38) of Lemma 1. The Hierarchy Theorem is proved. 

Assuming that w(0) is finite, we have for any n 

0<£ n <w(0). (3.110) 
Therefore, each of the monotonic sequences 

E x < £ 2 < £3 < ■ ■ ■ , in (A) 
£1 < £ 3 < £ 5 < • • • , in (B) 

and 

£2 > £4 > £e > • • • , in (B) 

converges to a finite limit £. By following the discussions in Section 5 of Ref.[4], one can 
show that each of the corresponding monotonic sequences of f n (x) also converges to a finite 
limit f(x). The interchange of the limit n — > 00 and the integrations in (3.13A) completes 
the proof that in Case (A) the limits £ and f(x) satisfy 

roc poo 

f(x) = l-2 <p- 2 (y)dy <p 2 (z)(w(z)-£)f(z)dz. (3.11L4) 

Jx Jy 

As noted before, the convergence in Case (A) can hold for any large but finite w(x), provided 
that w'(x) is negative for x > 0. In Case (B), a large w(x) may yield a negative f n (x), in 
violation of (3.23B) . Therefore, the convergence does depend on the smallness of w(x). One 
has to follow discussions similar to those given in Ref.[3] to ensure that the limits £ and f(x) 
satisfy 

f(x) = 1-2 f ' $-\y)dy f V <f ) 2 (z)(w(z)-£)f(z)dz. (3.111S) 

JO JO 
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4. Asymmetric Quartic Double-well Problem 



The hierarchy theorem established in the previous section has two restrictions: (i) the limitation 
of half space x > and (ii) the requirement of a monotonically decreasing perturbative 
potential w(x). In this section, we shall remove these two restrictions. 

Consider the specific example of an asymmetric quadratic double-well potential 

V{x) = l -g\x*-iy + g\x (4.1) 

with the constant A > 0. The groundstate wave function tp(x) and energy E satisfy the 
Schroedinger equation 

(T + V(x))i/;(x) = Ei/>(x), (4.2) 

where T = — as before. In the following, we shall present our method in two steps: We 
first construct a trial function cf)(x) of the form 

, , . f <f>+(x) for x > Q , 
4>{x) = { [,_) { , (4.3) 
^ y ' y <p [x] for x < 0. v ' 

At x = 0, 4>(x) and <p'(x) are both continuous, given by 

0(O) = + (O)=0-(O) (4.4) 

and 

0'(O) = +/ (O) = <T'(0) = 0, (4.5) 

with prime denoting J^, as before. As we shall see, for x > 0, the trial function 4>(x) = <p + (x) 
satisfies 

(T + V(x) + v + (x))<f) + (x) = E+<p + {x) (4.6a) 

with 

v' + (x) < 0, (4.7a) 
whereas for x < 0, 4>(x) = 4>~(x) satisfies 

(T + V(x) + v4x))<j>-(x) = Eq<P-{x) (4.66) 

with 

v'_(x) > 0. (4.76) 
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Furthermore, at x = ±00 

v + (oo) = v-(-oo) = 0. (4.8) 

Starting separately from <p + (x) and (f)~(x) and applying the hierarchy theorem, as we shall 
show, we can construct from <f)(x) another trial function 

( \ — j X + ( x ) f° r x > (A Q x 

X[X) ~ \ x~(x) for x<0 l4 ' yj 

with x{ x ) ar| d x'( x ) both continuous at x = 0, given by 

x (0) = X + (0) = x -(0) (4.10) 

and 

X'(0) = x + '(0) = x"'(0) = 0. (4.11) 
In addition, they satisfy the following Schroedinger equations 

(T + V(x)) X + (x) = E + X + ( x ) for a; > (4.12) 

and 

(T + V(x))x~(x) = E~x~( x ) for x < 0. (4.13) 

From V(x) given by (4.1) with A positive, we see that at any x > 0, V(x) > V(—x); therefore, 
E + > E~. 

Our second step is to regard x( x ) as 3 new trial function, which satisfies 

(T + V(x)+w(x))x( x ) = E oX (x) (4.14) 
with w(x) being a step function, 

»w4<">'{" (4l5) 

and 

E = i(E + + E-). (4.16) 

We see that w(x) is now monotonic, with 

w'(x) < (4.17) 

for the entire range of x from —00 to +00. The hierarchy theorem can be applied again, and 
that will lead from x( x ) to ip(x), as we shall see. 
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4.1 Construction of the First Trial Function 

We consider first the positive x region. Following Sec. 2.1, we begin with the usual 
perturbative power series expansion for 

rj}{x) = e~ 9S{x) (4.18) 

with 

gS(x) = gS (+) + 5x(+) + g~ 1 S 2 (+) + ■■■ (4.19) 

and 

E = gE {+) + + (T^W + • • • , (4.20) 

in which S n (+) and E n (+) are ^-independent. Substituting (4.18)-(4.20) into the Schroedinger 
equation (4.2) and equating both sides, we find 

S>(+) = x'-l, (4.21) 

S'o(+)S[(+) = ±SZ(+) + Xc-E (+), (4.22) 

etc. Thus, (4.21) leads to 

5o(+) = i(x-l) 2 (x + 2). (4.23) 

Since the left side of (4.22) vanishes at x — 1, so is the right side; hence, we determine 

E (+) = 1 + A, (4.24) 

which leads to 

5i(+) = (l + A)ln(l + x). (4.25) 

Of course, the power series expansion (4.19) and (4.20) are both divergent. However, if we 
retain the first two terms in (4.19), the function 

/ 1 \ 1+A 

0(+) = e-gSoM-SiW = e' 9 ^ (4.26) 

serves as a reasonable approximation of tp(x) for x > 0, except when x is near zero. By 
differentiating (/>(+), we find 0(+) satisfies 

(T + V(x) + u + (x))<f>(+) = 9 E (+)<f>(+) (4.27) 

where 

«,(») = | (g,(+)'-af(+)) = (i y^ A) - («8) 
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In order to construct the trial function <f>(x) that satisfies the boundary condition (4.5), 
we introduce for x > 0, 



}(+)= ( T l_y +A e -|^(+), (4.29) 



and 



, +f , 1 I (g + l + \)(f>(+) + (g-l-\)<f>(+), for0<x<l 

{X) = Tg ■ \ ((g + 1 + A) + (g - 1 - A)e"S«) 0(+), for x > 1 (430) 

so that + (x) and its derivative +/ (x) are both continuous at x — 1, and in addition, at 
x = 0we have +/ (O) = 0. For x < 0, we observe that V(x) is invariant under 

A -> -A and a; -> -x. (4.31) 

The same transformation converts + (x) for x positive to 4>~(x) for x negative. Define 

,- ( v 1 / (ff + 1 - A)0(-) + - 1 + A)^(-), for -Kx<0 
(x) = 2^ " \ ((«? + 1 - A) + (g - 1 + \)e-b) *(-), for x < -1 (432) 



where 



0(-) = (t^) 1 V**" (4-33) 

0(-) = {j^j A e-3^ sSo( - ) (4.34) 

and 

5 (-) = i(x + l) 2 (-x + 2). (4.35) 

Both _ (a;) and its derivative <f)~'(x) are continuous at x = —1; furthermore, + (x) and (f)~(x) 
also satisfy the continuity condition (4.4) and (4.5), as well as the Schroedinger equation (4.6a) 
and (4.6b), with the perturbative potentials v+(x) and i>_(x) given by 

v + (x) = u + (x) + g+(x) (4.36a) 

and 

v_ (x) = m_ (x) + g_ (x) , (4.366) 
in which u + (x) is given by (4.28), 
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. , Mszl=m±^M s^fgl , for < x < 1 

0+0*0 = { (g+l+A)+(g-l-A) e -TJ«+ 2 ^o(+) (4.38a) 

, for i > 1 

and 

f 2 fl (g-l+A)(l-A-As) e-l^+W^ , _ 

£-(x)= ( 9+ l-A) + ( 9 -l+A) e -§»+^0(-) ' KX<U (43g&) 

{ . for x < -1 

In order that <?+(x) be positive for x > and U-(x), g~(x) positive for x < 0, we 
impose 

A<1 and 1 + A < g, (4.39) 
in addition to the earlier condition A > 0. From (4.28) and (4.37), we have 

V( x ) = - (1 + A)(2 + A ) <0 for x>0 (4.40a) 

+ v ' (1 + x) 3 v 7 

and 

u'_(x)= ( X + A )( 2 + A ) >0 for x<0 . (4.406) 

(1 — x) d 

Likewise, from (4.38a) and (4.38b), we find 

~> , ^ f (-JTTT-r 2 g (i-x^)( fl +i+A) ) g ( x ) f or 0<x<l ,. . n , 

5f+(x) = < V 1+A - Ax (p + i+A)+( ff -i-A) e - 3 9+ 2 » s o(+)y »+v y> (4.41a) 
[ , for x > 1 



and 



-^ + — , 2g ! 1 ": 2Kg+1 ^L Sn( -J ^(^), for -Kx<0 



§L(x) = \ \ i— A— Ax -r (9+1 _ A)+(9 _ 1+A) e -f 9+ 2 9So( -); - x^^^w (4 416) 

[ , for x< -1. 

Furthermore, as x — > ±1, 

saw - - (g+1 2 : ( l ) ; 1 ( ;- > 1 :t) e -i/ '- 1 ' « 4 - 42a > 

and 

lim (x) = 2g(g ~ 1 + A) 6 " |g -S(x + 1). (4.426) 

— - 1 (g + l-X) + {g-l + X) e-h 

Thus, for x > 0, we have 

g' + (x) < (4.43) 
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and, together with (4.36a) and (4.40a), 

v' + (x) < (4.44a) 
for x positive. On the other hand for x < 0, g'_(x) is not always positive; e.g., at x — 0, 

#-(0) = + (9 + 1 - A)) (1 - X)(g - 1 + A) 

which is positive for g > 3A ^~ a2 , but at x = — 1+, 



\9 



#_(-l+) = -A£_(-l+ = ^ ^ < 0. 

(g + l-\) + (g-l + \)e-i» 

However, at x — — 1, u'_(— 1) = |(1 + A)(2 + A). It is easy to see that the sum u_ + g_ = i>_ 
can satisfy for rr < 0, 

v'_(x) > if # » 1. (4.446) 

To summarize: + (x) and (f)^(x) satisfy the Schroedinger equation (4.6a) and (4.6b), with 
v±(x) given by (4.36a) and (4.36b), 

E+ = gE (+) = g(l + A) (4.45a) 

and 

E~ = gE (-) = g(l - A) (4.456) 
and the boundary conditions (4.4) and (4.5). In addition, v±(x) satisfies 

v+ (oo)=0, v _(-oo) = (4.46) 
and the monotonicity conditions (4.7a) and (4.7b). 
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4.2 Construction of the Second Trial Function 

To construct the second trial function X (x) introduced in (4.9), we define f ± (x) by 

x +(x) = (f) + (x)f + (x) for x > (4.47a) 

and 

x -(x) = <f>-(x)f-(x) for x < (4.476) 

To retain flexibility it is convenient to impose only the boundary condition (4.11) first, but 
not (4.10); i.e., at x = 

X +/ (0) = X-'(0) = 0, (4.48) 

but leaving the choice of the overall normalization of x + (0) and x~(0) t0 be decided later. 
We rewrite the Schroedinger equations (4.12) and (4.13) in their equivalent forms 

(T + V(x) + v+(x) - E^)x + (x) = (v+(x) - E + )x + (x) (4.49a) 

and 

(T + V(x) + v-(x) - E )x~{x) = (v-(x) - £~)x~{x) (4.496) 

where 

£ + = E£ - E + (4.50a) 

and 

£- = Eq-E-. (4.506) 

Because at x — 0, + '(O) = <f>~'(0) = 0, in accordance with (4.5), we have, on account of 
(4.48), 

/ + '(0) = /-'(0) = 0. (4.51) 
So far, the overall normalization of f + (x) and f~(x) are still free. We may choose 

/ + (oo) = 1 and /"(-oo) = 1. (4.52) 

From (4.6a), (4.47a), (4.49a) and (4.50a), we see that f + (x) satisfies the integral equation 
(for x > 0) 

oo oo 

f + {x) = 1-2 J ^ + {y)Y 2 dy J {^{z)f{v + {z) - £ + )f + (z)dz. (4.53a) 
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Furthermore, from (4.6a) and (4.49a), we also have 

oo 

J {<f) + {x)) 2 {v + {x) - £ + )f+(x)dx. (4.54a) 
o 

Likewise, f~(x) satisfies (for x < 0) 

x y 

f-(x) = 1 - 2 J (r(y)r 2 dy J (0») 2 M^) - £~)f-(z)dz (4.536) 

— oo — oo 

and 

o 

J ((f)-(x)) 2 (v-(x) - S~)r(x)dx. (4.546) 

— oo 

The function f + (x) and f~(x) will be solved through the iterative process described in Sec- 
tion 1. We introduce the sequences {fn( x )} ar| d } f° r ^ = 1, 2, 3, • • •, with 

oo oo 

= 1-2 J (<P + (y))- 2 dy J (<p + (z)) 2 (v + (z) - S^f^dz. (4.55a) 

for x > 0, and 

/"(*) = 1 " 2 / (riyT'dy J {<P~{z)) 2 {v^z) - E-)f~_ x {z)dz (4.556) 

— oo — oc 

for x < 0, where £^ satisfies 

oo 

J (<P + (x)) 2 (v + (x) - S^ft^dx (4.56a) 
o 

and 

o 

j {<t>-{x)) 2 {v-{x) - S-)f-^(x)dx. (4.566) 

— oo 

Thus, (4.55a) and (4.55b) can also be written in their equivalent expressions 

x y 

= /„ + (0) - 2 / (<f> + (y))- 2 dy J (0 + (z)) 2 K(z) - #)/» + -i(*)d*. (4.57a) 



for x > 0, and 



/„» = /n (0) - 2 / (0-(y))- 2 dy / (<T (*)) 2 (M*) - £-)f-^(z)dz (4.576) 



x y 



x y 

for a; < 0. 
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For n = 0, we set 

f+(x) = / "(x) = 1; (4.58) 

through induction and by using (4.55a)-(4.56b), we derive all subsequent f„(x) and 
Because f±(x) satisfies (4.44a), (4.44b) and (4.46), the Hierarchy theorem proved in Section 3 
applies. The boundary conditions / + (oo) = /~(— oo) = 1, given by (4.52), lead to /+(oo) = 
f~(-oo) = 1, in agreement with (4.55a) and (4.55b). According to (3.24)-(3.27) of Case (A) 
of the theorem, we have 

St < £ 2 + < £ 3 + < ■ ■ ■ , (4.59a) 
£f < £ 2 _ < £ 3 _ < • • • , (4.596) 
1 < f+(x) < f}{x) < f+(x) <■■■ (4.60a) 

at all finite and positive x, and 

1< fr(x) < / 2 -(s) < / s -(x) < • • • (4.606) 
at all finite and negative x. Since 

5+ < u+(0) (4.61a) 

and 

5-<u_(0) (4-616) 

with both f±(0) finite, 

lim £+ = £+ and lim S~ = S~ (4.62) 

both exist. Furthermore, by using the integral equations (4.55a)-(4.55b) for f„(x) and by 
following the arguments similar to those given in Section 5 of Ref. 13, we can show that 

Urn /+(*) =/+(*) and Bm /"(*) = /"(*) (4.63) 

also exist. This leads us from the first trial function 4>(x) given by (4.3) to 4>f + (x) and 4>f~(x) 
which are solutions of 

(T + V(x) -E + )(pf + (x) =0 for x > (4.64a) 

and 

(T + V(x)-E-)(f>f-(x)=Q for x < (4.646) 

with 

E + = E+ -S + and E~ = Eq - £" , (4.65) 
and the boundary conditions at x — 0, 

0'(O) = /+'(0) = /-'(0) = 0. (4.66) 

An additional normalization factor multiplying, say, f~(x) would enable us to construct the 
second trial function x(x) that satisfies (4.9)-(4.13). 
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4.3 Symmetric vs Asymmetric Potential 

As we shall discuss, the general description leading from the trial function x( x ) t0 the 
final wave function ip(x) that satisfies the Schroedinger equation (4.2) may be set in a more 
general framework. Decompose any potential V(x) into two parts 

- { ^1; £ ;io° 

Next, extend the functions V a (x) and Vfe(x) by defining 

V a (x) = V a {-x) for x < 

and (4.68) 

Vft(a;) = V^-ar) for x > 0. 

Thus, both V a (x) and Vb(x) are symmetric potential covering the entire rr-axis. Let Xa(x) 
and Xb(x) be the groundstate wave functions of the Hamiltonians T + V a and T + V b . 

(T + K(z))Xa(z) = £ a x«(^) (4.69a) 

and 

(T + V b (x)) Xb (x) = E bXb (x). (4.696) 
The symmetry (4.68) implies that 

Xa(x) = Xa(-x), Xb(x) = Xb(-x) (4.70) 

and at x = 

X'a(0) = X' 6 (0) = 0. (4.71) 

Choose the relative normalization factors of Xa and Xb, so that at rr = 

Xa(0) = Xb(0). (4.72) 

The same trial function (4.9) for the specific quartic potential (4.1) is a special example of 



with 



Xa{x) = X + {x) for x > 

and (4.74) 

X&(^) = X~(^) for x < 0. 
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In general, from (4.69a)-(4.69b) we see that x{ x ) satisfies 

(T + V(x) + w(x)) X (x) = E oX {x). (4.75) 
Depending on the relative magnitude of E a and E b , we define, in the case of E a > E b 

and 

E = E a ; (4.77a) 

otherwise, if E b > E a , we set 

A / x \ E b — E a for rr > ,. „ n 

W & = {o for x < ^ 

and 

= E b . (4.776) 

Thus, we have either 

w(oo) = and w'(x) < (4.78a) 

at all finite x, or 

w(-oo) = and w'(x) > (4.786) 

at all finite x. A comparison between (4.9)-(4.17) and (4.73)-(4.77a) shows that w(x) of 
(4.14) and the above w(x) differs only by a constant. 

As in (4.2), ip(x) is the groundstate wave function that satisfies 

(T + V(x))ij(x) = Eij(x), (4.79) 

which can also be written in the same form as (1.14) 

(T + V(x) + w(x) - E )ip{x) = (w(x) - £)ip{x), (4.80) 

with 

E = E -S. (4.81) 

Here, unlike (1.32), V(x) can now also be asymmetric. Taking the difference between i^(x) 
times (4.75) and x( x ) times (4.80), we derive 



/CO 
X(x)if>(x) (w(x) - S)dx = 0. (4.82) 
-co 



41 



Introduce 

= X{x)f(x), (4.83) 

in which f(x) satisfies 

roc yco 

f(x) = /(oo) - 2 / X -\y)dy X 2 (z)(w(z)-S)f(z)dz. (4.84) 
On account of (4.82)-(4.83), the same equation can also be written as 

f(x) = /(-oo) - 2 f X - 2 (y)dy f X 2 (z)(w(z) - S)f(z)dz. (4.85) 

J — oo J —oo 

Eq. (4.80) will again be solved iteratively by introducing 

Mx) = X(x)fn(x) (4.86) 
with tp n and its associated energy £ n determined by 

(T + V{x) + w(x) - Eo)M%) = («Kz) - 4)^-i(z) (4.87) 

and 

/oo 

X (x)Vn-iW(^^) - £„)<*r = 0. (4.88) 
-oo 

In terms of f n (x), we have 

/•oo roo 

f n (x) = f n (oo) - 2 / ;T%)*/ / x 2 (z)(w(z) - £ n )f n ^(z)dz. (4.89) 
On account of (4.88), we also have 

/oo 

X 2 (rr)(wOr) - £ n )fn(x)dx = (4.90) 
-oo 

and 

fn(x) = /n(-oo) - 2 f x~ 2 (y)^ f X 2 (z)(w(z) - £ n )f n ^{z)dz. (4.91) 

J —oo J —CO 

For definiteness, let us assume that 

E a > E b (4.92) 

in (4.69a)-(4.69b); therefore w'(x) < and w(oo) = 0, in accordance with (4.76a). Start 
with, for n — 0, 

/ o (a0 = 1, (4.93) 
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we can derive {E n } and {f n (x)}, with 

E n = E - £ n , (4.94) 
by using the boundary conditions, either 

/„(oo) = l for all n, (A) 

or 

/ n (-oo) = 1 for all n. (B) 

It is straightforward to generalize the Hierarchy theorem to the present case. As in Section 3, 
in Case (A), the validity of the Hierarchy theorem imposes no condition on the magnitude of 
w(x). But in Case (B) we assume w(x) to be not too large so that (4.91) and the boundary 
condition /„(— oo) = 1 is consistent with 

fn(x) > (4.95) 

for all finite x. From the Hierarchy theorem, we find in Case (A) 

E 1 > E 2 > E 3 > ■ ■ ■ (4.96) 



and 



while in Case (B) 



and 



I < fi(x) < f 2 (x) < f 3 (x) < ■ ■ ■ , (4.97) 



E x > E 3 > E 5 > ■ ■ ■ , (4.98) 
E 2 < E 4 < E 6 < ■■■ , (4.99) 
l</i(z)</ 3 (x)</ 5 (a;)<-.. (4.100) 



1 > f 2 (x) > h{x) > f 6 (x) > ■ ■ ■ . (4.101) 



A soluble model of an asymmetric square-well potential is given in Appendix A to illustrate 
these properties. 
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5. The N- Dimensional Problem 

The iV-dimensional case will be discussed in this Section. We begin with the electrostatic 
analog introduced in Section 1. Suppose that the (n — l)th iterative solution / n _i(q) is 
already known. The nth order charge density o"„(q) is 

cr„(q) = (w(q) - £ n )0 2 (q)/„_i(q), (5.1) 

in accordance with (1.23)-(1.24). Likewise, from (1.26) and (1.29) the dielectric-constant k 
of the medium is related to the trial function 0(q) by 

><q) = 2 (q), (5.2) 

and the nth order energy shift S n is determined by 

J a n (q)d N q = 0. (5.3) 

In the following we assume the range of w(q) to be finite, with 

w(oo) = (5.4) 

and 

< w(q) < W rnax . (5.5) 

Introduce 

K(W) = J K(q)5(w(q) - W)d N q, (5.6) 

where 5(w(q) — W) is Dirac's 5-function, IT" is a constant parameter and the integrations in 
(5.3) and (5.6) are over all q-space. Similarly, for any function F(q), we define 

F av (W) = KiW)- 1 J F(q)K(q)5(w(q) - W)d N q. (5.7) 

In the iV-dimensional case, the generalization of [F], introduced by (3.15), is 

[F] = J F(q)0 2 (q)^g. (5.8) 

In terms of F av (W), (5.8) can also be written as 

[F] = / F av (W)K(W)dW. (5.9) 
Jo 

Thus from (5.1) and (5.3) we have 

£n = Kq)/n-i(q)]/[/n-i(q)], (5.10) 
the n-dimensional extension of (3.14). 
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Following (1.27)-(1.28), the nth order electric field is — |V/ n and the displacement field 

is 

D n = ~kVf n . (5.11) 
The corresponding Maxwell equation is 

V • D n = a n . (5.12) 

Eqs.(5.11) and (5.12) determine /„ except for an additive constant, which can be chosen by 
requiring 

minimum of / n (q) = 1. (5.13) 

Therefore, 

/n(q) > 1. (5.14) 

As in the one dimensional case discussed in Section 3, (5.10) gives the same condition of fine 
energy tuning at each order of iteration. It is this condition that leads to convergent iterative 
solutions derived in Section 3. 
We now conjecture that 

lim E n = E = E - E (5.15) 

n^oo 

and 

lim / n (q) = /(q) = ^(q)/0(q). (5.16) 

also hold in higher dimensions. Although we are not able to establish this conjecture, in the 
following we present the proofs of the iV-dimensional generalizations of some of the lemmas 
proved in Section 3. 

Lemma 1. For any pair / m (q) and //(q) if at all W within the range (5.5), 



and 

d (f%{wy 



dW \ ft v (W)j 
Proof . For any function !F(W), define 



> 0, then £ m+1 > £ l+1 , (5.17) 



< 0, then £ m+1 < £ l+1 . (5.18) 



ax 

< F(W) >= J K(W)F(W)dW, (5.19) 
o 

Thus for any function F(q), we have 

[F(q)] =< F av (W) >; (5.20) 

therefore, 

E n [fn-M)]=£n<f: V - 1 {W)> (5.21) 
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and 

Kq)/n-i(q)] =< WfZ.iW) > . (5.22) 
By setting the subscript n in (5.10) to be m + 1, we obtain 

£ m+ i < >=< TO > • ( 5 -23) 

Also by definition (5.19), 

S+i < >=< > • (5-24) 
The difference of (5.23) and (5.24) gives 

(£ m+1 - S l+1 ) < f™{W) >=< (W - £i + i)f™{W) > . (5.25) 
From (5.10) and setting the subscript n to be I + 1, we have 

0=<(W-S l+1 )fr(W)>. (5.26) 

Regard f? v {£i+i) and /™(£ J+ i) as two constant parameters. Multiply (5.25) by f? v (£i +1 ), 
(5.26) by f™(£ i+1 ) and take their difference. The result is 

(£ m+ i-£i + i)fr(£i + i) < m >=< (^-^iK/rfei/rw-zrw/:^) >, 

(5.27) 

analogous to (3.43). 

(i) "f iK/r W/ifW) > 0, then for W < £ m 

Thus, the function inside the bracket < > in (5.21) is positive, being the product of 
two negative factors, (W - £ l+1 ) and {f? v {£ l+1 )f™{W) - fr(W)f™(£ l+1 )). Also, when 
W > £i+±, these two factors both reverse their signs. Consequently (5.17) holds. 

00 lf mifm(W)/fr(W)) < °- we see that for W < ( 5 - 28 ) reverses its sign, and 

therefore the function inside the bracket < > in (5.27) is now negative. The same negative 
sign can be readily established for W > £ t+1 . Consequently, (5.18) holds and Lemma 1 is 
established. 

Lemma 2. Identical to Lemma 2 of Section 3. 

In order to establish the iV-dimensional generalization of Lemma 3 of Section 3, we define 

Q n (W) = J dW 1 J a n (q)5(w(q) - W 1 )d N q. (5.29) 
w 

Because of (5.3), Q n (W) is also given by 

w 

Qn(W) = -JdW 1 J a n (q)5(w(q) - W x )d N q. (5.30) 
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We may picture that the entire g-space is divided into two regions 

I. w(q) > W (5.31) 

and 

//. w(q) < W, (5.32) 

with Q n (W) the total charge in /, which is also the negative of the total charge in II. By 
using (5.1) and (5.7), we see that 

~^£P- = ( W - Wi(W»0- (5-33) 
Lemma 3. For any pair / m (q) and /;(q) if at all W within the range (5.5) 

,,x d (gj)\ d ( Q m+1 (W) \ 

w dw[jrW)) >0 then M^wJ >0 (5 ' 34) 
(,,) dw{fr(w)J dw{ Ql+1 (w)) <0 - (5 - 35) 

Proof Note that (5.34) and (5.35) are very similar to (3.56) and (3.57). As in (3.60), define 

Z = Qi+i(W) and v = Q m+1 (W). (5.36) 



From (5.33), we have 



and 



^'""'"^^-U/rW^) (5-37) 



dW 
dQi +1 (W) 



dW 



(W-£ l+1 )fr(W)K(W). (5.38) 



Therefore, 



where 



Furthermore, 



dr, _ dQ m+1 (W) , dQ l+1 (W) _ ( r o Q] 

dt ~ dW / dw ~ { ] lrW) { ] 

r iW) = ^|±i . (5.40) 



d fdr]\ dQi +1 \ d ( dQ m+1 

di \d£J \ dW ) dW\dw'dW) \ dW ) dW \ ' ff v 



'dQi+iY 1 (dr_f^,d (f%\\ 

^ dw ) \dw fr dw \ fr ) )' 1 ■ ' 
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where 

dr _ £ m +i - Si + i (kao\ 

dW~(W-£ l+i r { ] 

analogous to (3.61)-(3.64). 

According to (5.30), at W = 

Qm+i(0) = Q J+ i(0) = (5.43) 

and according to (5.29), at W = W max 

Q m+ i(W max ) = Q l+1 (W max ) = 0. (5.44) 

From (5.37), we see that the derivative -^Q m+ i(W) is positive when W < £ m +i, zero at 
W = £ m +i and negative when W > £ m +i- Likewise, from (5.38), -^Qi + i{W) is positive 

when W < £i+\, zero at W = £i+\ and negative when W > £i +1 . Their ratio determines ^ 

(i) "f Mf^/fi V ) > °. from lemma L we have 

£ m +i > £i+i (5.45) 

and therefore, on account of (5.42) 

(IV 

" >0. (5.46) 



dW 
At W = 0, 

r(0) = §^ > 1. (5.47) 

As V4 7 increases, so does r(W). At W — £\+\, r(W) has a discontinuity, with 

r(£ I+1 -) = oo (5.48) 

and 

r(£ l+1 +) = -oo. (5.49) 
As W increases from £i + i, r(W) continues to increase, with 

r(£ m+1 ) = (5.50) 

and 

r(W max ) = W ™ x - S ? +1 < 1. (5.51) 
It is convenient to divide the range < W < W max into three regions: 

A. < W < £ l+1 (5.52) 

B. £ l+1 <W< £ m+1 (5.53) 
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C. S m+1 < W < W max . (5.54) 

Assuming -^{fm I fi v ) > ®< we sna " snow separately -^(Q m+1 /Qi +1 ) > in these three 
regions. 

In region B, Qz+i is decreasing, but Q m+ i is increasing. Clearly, 

' (Qm+i/Qi+i) > 0. (5.55) 



In region A, ^§±1 > q, r(W) is positive according to (5.47)-(5.48) and is always > 
from (5.46). Therefore from (5.41), 



-j± > in A. (5.56) 



In region C, < 0, but r(W) and are both positive. Hence, 

/7 <0 in C. (5.57) 



Within each region, rj = Q rn +i{W) and £ = Q;+i(W / ) are both monotonic in V4 7 ; therefore 
?7 is a single-valued function of £ and we can apply Lemma 2 of Section 3. 

In region A, at W = both Q m +i(0) and Q;+i(0) are according to (5.43), but their 
ratio is given by 

Qm+i(0) _ f dQ m+1 (W) , dQ l+1 {W) \ 

Q^W" " I dw I -^^) w ^ ■ (5 - 58) 

Therefore 

( ^£ ) vf=o ( £ ) w=o 

Furthermore, from (5.56), ^ > 0. It follows from Lemma 2 of Section 3, the ratio rj/^ is an 
increasing function of £. Since 



dW dW 

as also have 



d (Q m+l \ di d M >0 in A (561) 



In region C, at W = W max , both Q m+ i(W max ) and <3/+i(W ma s) are according to (5.44). 
Their ratio is 

Q m +i(W max ) _ ( dQ m+l (W) j dQi +1 {W) \ . 2 

Qi + i(w max ) { dw I dw ) w= ' 1 • } 
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Which gives at W = W„ 



?7 



(5.63) 



As W decreases from W„ 



to £ m +i in region C, since £ m+ i > £i + i, we have 
d£ dQi +1 



dW dW 



< in C. 



(5.64) 



Furthermore, from (5.57), |? < in region C. It follows from Lemma 2 of Section 3, the 



ratio 77/^ is a decreasing function of £, which together with (5.64) lead to 

d I Qm+l \ 



dW \ Q l+1 J 



> in C. 



(5.65) 



Thus, we prove case (i) of Lemma 3. Case (ii) of Lemma 3 follows from case (i) by an 
exchange of the subscripts m and I. Lemma 3 is then proved. 

So far, the above Lemmas 1 and 3 are almost identical copies of Lemmas 1 and 3 of 
Section 3, but now applicable to the iV-dimensional problem. Difficulty arises when we try to 
generalize Lemma 4 of Section 3. 

It is convenient to transform the Cartesian coordinates qi, q 2 , ■ ■ ■ , qN to a new set of 
orthogonal coordinates: 



with 
and 

where i or j — 1, 2 



(91, Q2, qn) -> Mq), A(q), ^jv-i(q)) 
Vw • va = 



VA • VPj = for ^ j, 
N — 1. Introducing 

hi = l/(Vw) 2 , h\ = 1/(VA) 2 

w = h w Vw and A = fojVA 
In terms of the new coordinates, the components of D n are 

(D n ) w = w-D n and (D n )i = fa ■ D n . 

Its divergence is 

f d N—i q 



(5.66) 

(5.67) 
(5.68) 

(5.69) 
(5.70) 

(5.71) 
(5.72) 
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W N-l 



Combining (5.12) with (5.30), we have 

w 

QniW) = - J dw J h w h {V ■ D n ) n dPi ; (5.73) 
o 1 

therefore, 

- JV-l 

QniW) = -f (D n ) w h p J] Wi, ( 5 - 74 ) 
J 1 

in which the integration is along the surface 

w(q) = W. (5.75) 

From (5.11) and (5.71), it follows that 

(D n ) w = = ~nh w Vw ■ V/„. (5.76) 

In terms of curvilinear coordinates, (5.7) can be written as 

. N-l 

F av (W) = K(W)~ l f F(q) K (q)h w hp JJ d & ■ (5-77) 
J i 

Substituting (5.76) into (5.74), we find 

Q n (W) = l -K{W){Vw ■ Vf n ) av . (5.78) 
Because h~ l (df n / dw) = w ■ V/ n = h w (Vw ■ V/ n ), (5.78) can also be written as 



av 



Qn(W) = 1 -K(W)( [ h~J^j . (5.79) 

Here comes the difficulty. While the above Lemma 3 transfers relations between I f\ v to 
those between Q m +i/Qi+i, the latter is 



av 



[ h « 8 -W) i{^ s -t) < 5 - 8o > 

which is quite different from ■£yfm+i/'£vft+i- This particular generalization of the lemmas 
in higher dimensions fails to establish the Hierarchy Theorem. 

For the one-dimensional case discussed in Section 3, we have w' < and x > 0; con- 
sequently (5.80) is fm+i/fl+i- Therefore, Lemma 4 of Section 3 can also be established by 
using (5.80), and the proof of the Hierarchy Theorem can be completed. 
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Appendix 



A.l A Soluble Example 

In this Appendix, we consider a soluble model in which the potential V(x) of (4.67) is 



V{x) = < 



oo 

\W 2 


oo 



for 



7 < x 

a < x < 7 
-a < x < a 
-7 < x < —a 
x < -7, 



with W 2 > fj 2 and 



7 = a + j3. 

Following (4.68), we introduce two symmetric potentials: 

V a (x) = V a (-x) and V b (x) = V b (-x) 



with, for x > 0, 

V a (x) 

and 

V b (x) 



00 

\w 2 



00 


\w 2 



for 



for 



7 < x 
a < x < 7 
< x < a 



7 < x 
a < x < 7 
< x < a, 



so that (A.l) can also be written as 

' V a (x) 



V(x) = 



V b (x) 



for 
for 



x > 
x < 0. 



Let ip(x), Xa(x) and Xb{x) be respectively the groundstate wave functions of 

(T + V(x))ip(x) = EiP(x) 

(T + V a (x)) X a(x) = E aXa (x) 

and 

(T + V b (x)) Xb (x) = E bXb (x). 
For \x\ > 7, since V(x) = 00, we have 

^0*0 = Xa(x) = Xft(^) = 0. 



(Al) 



(A2) 
(A3) 

(A4) 



(A5) 



(A6) 

(A7) 
(A8) 

(A9) 
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For |x| < 7, these wave functions are of the form 



ip(x) oc < 



sin k (—x + 7) 
cosh q(x — 5) 
ship (x + 7) 



a < x < 7 
for —a < x < a 

— 7 < x < —a, 



Xa(x) OC < 



sin k a (-x + 7) 
coshg a x 
sin k a (x + 7) 



for 



and 



Xb{x) oc < 



sinp 6 (-x + 7) 
coshg b x 
k sinp 6 (a; + 7) 



for 



a < x < 7 
—a < rr < a 
—7 < rr < —a 



a < x < 7 
-a < x < a 
-7 < rr < —a. 



By substituting these solutions to the Schroedinger equations (A.7)-(A.9), we derive 

E = \p 2 = \(W*- q *) = \^ + k% 



E n 



and 



\{w"- q l) = \{^ + kl) 



1 j 1 



The continuity of ■?///■?/> at x = ±a relates 

— k(5 cot k(5 = q/3 tanh q (a — 5) 

and 

— pf3 cot p/3 = q(3 tanh 5 (a + 5) 

with /3 given by (A. 2). In the following, we assume the barrier heights W and 

W = (W 2 - /1 2 ) 1 / 2 



(A10) 



(All) 



(A12) 



(A13) 
(A 14) 

(A15) 

(A16) 
(A17) 
(A18) 



to be much larger than k and p; therefore, the wave function ip is mostly contained within the 
two square-wells; i.e., k(3 and p(3 are both near n. We write 



kp = n-9, pp = n-9 

and expect 9 and 9 to be small. Likewise, introduce 

k a (3 = n - 9 a , p b p = 7r - 9 b . 



(A19) 



(A20) 
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The explicit forms of these angles can be most conveniently derived by recognizing the 
separate actions of two related small parameters: one proportional to the inverse of the barrier 
height 

(Wfi)- 1 « 1 (A21) 

and the other 

e -2Wa <<< 1 ( A22 ) 

denoting the much smaller tunnelling coefficient. 

To illustrate how these two effects can be separated, let us consider first the determination 
of 9b given by (A. 20). The continuity of x'b/Xb at x — ±ct gives 

— Pbj3 cot pb(3 = qbP tanh q b a. (A 23) 

From (A. 15), we also have 

W 2 =pl + ql (A24) 

Although the two small parameters(A.21) and (A. 22) are not independent, their effects can 
be separated by introducing and that satisfy 

- poo/3 cot poo/3 = ( A25) 

and 

W 2 = pi + (A26) 

Physically, p^ and q^ are the limiting values of p b and q b when the distance 2a between the 
two wells — > oo, but keeping the shapes of the two wells unchanged. Hence (A. 23) becomes 
(A.25). Let 

Poo/3 = 7T - #oo • (A27) 

From (A.25), we may expand 9^ in terms of successive powers of (W/3)^ 1 : 

9 ^wff( 1 -m +{1+ ^ )( m ?+0( m f ) (A28) 

which determines both poo and goo- By substituting 

9 b = 9 oa + v x e~ 2q °° a + 0{e- 4q °° a ) (A29) 
into (A. 23) and using (A.24)-(A.28), we determine 

(Pooqoo\ 2goo/3 

= ^JTT- <AM) 
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Likewise, the continuity of x' a /Xa at x — ±a gives 

— k a (5 cot k a (5 = q a f3 tanh k a a, (A 31) 

with 

W 2 = W 2 -fx 2 = k 2 a + q 2 a . (A32) 
As in (A. 25), we introduce koo and that satisfy 

- kvo/3 cot /coo/9 = goo/3 (A. 33) 

and 

W 2 = kl + ql. (A34) 
Similar to (A.27)-(A.28), we define 

kooP = 7T - ^oo . (A35) 

and derive 

§ 0O = ^-(l-J- + (l + —)(J-) 2 + O(i) 3 ^ . (A36) 

As in (A.29)-(A.30), we find # a = n - k a (3 to be given by 

9 a = 9 00 + z>x e~ 2 ^ a + 0(e~ 4 ^ a ) (A37) 

with 

= / ^ogoo\ J|^L ( } 

V W 2 J qooP + 1 V 7 

To derive similar expressions for 9 and of (A. 19), we first note that the transformation 

+ 5 (A39) 
brings (A. 23) to (A. 17), provided that we also change 

qb -> q, Pb^P 

and therefore 

9 b -> 0. (A40) 
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Since according to (A.l), the asymmetry of V(x) is due to the term |yU 2 > in the positive 
x region, it is easy to see that 

5 > 0, (A41) 

as will also be shown explicitly below. Thus, from (A. 29) and through the transformations 
(A.39)-(A.40), we derive 

9 = #oo + Oi (A42) 

where 

e 1 = V\ e - 2 «~( Q+<5 ) + 0{e-^ {a+6) ) (A43) 
with v\ given by (A. 30). Likewise, we note that the transformation 

a -> a - 5 (AAA) 
brings (A. 31) to (A. 16), provided that we also change 

k a k, q a — > q 

and therefore 

6 a -> 6. (A45) 
Here, we must differentiate three different situations: 

(i) a > 5, 

(ii) a = 5 (A46) 

and 

(iii) a < 5. 

In case (i), when 

e -2 qoo (a-S) << 1 ( AA? j 

from (A. 37) and through the transformations given by (A.44)-(A.45), we find 

6 = L + h (A48) 

where 

0! = f>! e" 2 ^( a - 5 ) + 0(e- 4 ^ (a - 5) ) (A49) 
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with z>i given by (A. 38). According to (A. 13) and (A. 19), we have 

= (7r-e O0 -e l ) 2 -(n-e oo -e l ) 2 , (aso) 

which leads to 

H 2 (3 2 + (2tt -On- D(0oo - L) = -2(tt - 9^)9, + 2(vr - L)0i +0\- 9j. (AM) 
Since in accordance from (A. 28) and (A. 36), we find 

«--*-I(iF-F)-?(^-Wi) + -- (^)«^ (ab, 

and 

^ + ^ = ^(i 7 + i)+...«2 7 r. (A53) 

Thus, the left side of (A. 51) is dominated by its first term, fi 2 f3 2 . Since 9\ and 9 1 are 
exponentially small, we can neglect 9\ — 9\ in (A. 51). In addition, because 9^ and 9^ are 
much smaller than 2n, (A. 51) can be reduced to 

li 2 (3 2 = 2n(9 1 - B x ) = 2tt (z>! e - 2 *»(°-«) - v x e - 2q °° (a+5 ^ (A54) 

which gives the dependence of 5 on fx 2 . It is important to note that an exponentially small jj 2 
can produce a finite 5. For 5 < a, at x = 5 we have, in accordance with (A. 10) 

if>'(6) = 0, (A55) 

which gives the minimum of ip(x). The wave function ip(x) has two maxima, one for each 
potential well. 

In case (ii), a = 5 and (A. 16) gives cot kf3 = 0, and 9 takes on the critical value 9 C with 

9 C = \. (A56) 

In case (iii), 9 > |, k/3 < f and ^(x) has only one maximum. 
As in (4.73) and (4.75) we introduce x( x ) through 

Y ( x ) = i Xa(x), for 0<x< 7 

X[X) -\ X b(x), for - 7 <a;<0 {A ^ <} 

so that 

(T + V{x) + = E oX (x), . (A58) 
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in which, same as (4.76a)-(4.77a), 

^ x) = { E ° a -E b for r -7<^<0 (A.59) 

and 

E = E a (AM) 

with E a and E b given by (A. 14) and (A. 15). Since according to (A.4)-(A.5), V a (x) > V b (x), 
we have 

E a > E b ; (AM) 

therefore, 

w'(x) < 0. (AM) 
Write the Schroedinger equation (A. 7) in the form (4.80): 

(T + V(x) + w(x) - E )ip(x) = (w(x) - £)ip(x) (AM) 

with 

E = E -S = E a -S. (AM) 

As in (4.82), we have 

n 

I x{ x )i ) { x ){w(x) — £)dx = 0. (A. 65) 

In all subsequent equations, we restrict the rr-axis to 

\x\ < 7, (A66) 
and set vp(x), x( x ) positive. Define 

f(x)=^(x)/ X (x). (AM) 
We have, as in (4.84)-(4.85), 

f(x) = /( 7 ) -2 [\- 2 (y)dy r x 2 (z)(w(z)-S)f(z)dz, (AM) 

J x Jy 

or, on account of (A. 65), the equivalent form 

f{x) = /(- 7 ) " 2 f X - 2 (y)dy [ V X 2 {z)(w(z) - S)f(z)dz. (AM) 

J— y J— *y 
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The derivation of f(x) is given by 

f{x) = -2 x - 2 (x)h(x) (A70) 

where 

h(x) = f X 2 (z)(w(z) - S)f(z)dz (A71) 
J-i 

or equivalently, 

h(x) = - r x 2 (z)(w(z)-S)f(z)dz. (A72) 

Jx 

In order to satisfy (A. 65) and by using (A. 59), we see that 

(E a -E b )>S>0. (A73) 

Thus, (A.71)-(A.73) give 

(-7) = M7) = 0, 

(A74) 

h'(x) > for - 7 < x < 

and 

The positivity of h(x) gives 

f{x) < for \x\ < 7. (A75) 

A.2 A Two-level Model 

Before discussing the iterative solutions for fix) and S, it may be useful to first extract 
some essential features of the soluble square-well example. Let us first concentrate on Case 
(i) of (A. 46), with the parameters a and 5 satisfying 

e -2 goc a <<; e -2 qoc (a-S) <<c ^ ( A7g ) 

We shall also neglect iWa)' 1 or when compared to 1. Thus, from (A.27)-(A.28), 

we have 

Poo = ; (A77) 
in addition, from(A.30) and (A. 38) we find 



K-i) 




Ki) = o, 




h(x) 


> 


for 


\x\ < 7, 


h\x) 


> 


for 


- 7 < x < 


h\x) 


< 


for 


< x < 7. 



0, 9* Ul C* ?Poo ^ ^7T_ 
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From (A. 54), we have 



I, ,2 ~ K~ -2 qoo (a-5) ~ 2?I " 2 -2 qoo (a-S) 



On account of (A.15), (A.20), (A.27) and (A.29), 

E b = \vl = liPo, - U je-^) 2 

which, for 

E = -p 2 

gives 

It? 2tt 2 

rp ~ rp roc -2q x a ~ rp -2q x a 

On the other hand, from (A. 13), (A. 19) and (A.42)-(A.43), we see that 

F _ I 2 ~ 1/ _ V ± p -2qoo(a+S)\2 ~ Z7 _ -2 goo {a+6) 

tj ~ 2 P ~ 2 [Pco P ~ W[33 

Thus, under the condition (A. 76), we find 

-/i 2 >> Eoo — Eb >> Eoo — E. 

As we shall see, these inequalities can be understood in terms of a simple two- 
Introduce 

A = Eqo — E b . 

We note that from (A. 82), 

2tt 2 

and from (A.79) and (A.83), 

A I /i 2 e - 2fe<5 and E^ - E = \e~ 2q °° s . 
Consequently, the three small energy parameters in (A. 84) are related by 

^ 2 (E 00 -E)^X 2 , 

From e~ 2qocS « 1 and (A. 76) , we see that 

1 2 2tt 2 
A << -/T << 



2 r ~ ' W(3 3 
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in accordance with (A. 79) and (A. 84). To understand the role of the parameter A, we may start 
with the definition of Vb(x), given by (A. 5), keep the parameters (3 = 7 — a and \W 2 fixed, 
but let the spacing 2a between the two potential wells approach 00; in the limit 2a — > 00, 
we have E b — > E^. Thus, A = E^ — E b is the energy shift due to the tunneling between the 
two potential wells located at x < —a and x > a in V b (x). 

There is an alternative definition for A, which may further clarify its physical significance. 
According to (A. 3), V b (x) is even in x; therefore, its eigenstates are either even or odd in x. 
In (A. 9), Xb( x ) is the groundstate of T + V b (x), and therefore it has to be even in x. The 
corresponding first excited state Xod is odd in x; it satisfies 



(T + V b (x))xod{x) = E od Xod{x). 



We may define A by 



(A90) 
(A91) 

and regard (A. 85) and (A. 86) both as approximate expressions, as we shall see. 

Multiplying (A. 9) by Xod(x) and (A. 90) by Xb(%), then taking their difference we derive 



2A = E od — E, 



~2 (Xod( x )Xb(x) - X b ( x )Xod(x))' = (E od - E b )xod(x)Xb(x). 



(A92) 



From (A. 12), we may choose the normalization of Xb so tnat 



\ (^)^P6(-x + 7) 



Xb(x) = Xb(-x) = < 



Correspondingly, 



Xod{x) = -Xod{-x) 



with 



cosh q b x 



for 
for 



sinh q od x 



for 
for 



a < x < 7 
< x < a. 

a < x < 7 
< x < a 



(A93) 



(A94) 



Eb = \vl 



and 



E od = -Pod- 



As in (A. 25) and (A. 26), q od and p od are determined by 

-PodP cot p od (3 = q od (3 coth q od a 

and 

W 2 = p 2 od + q 2 od . 



(A95) 

(A96) 
(A97) 
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At x = 0, we have 

Xod(0) = 0, X b(0) = 1 

XoM = q od and ^(0) = 0. 

Integrating (A. 92) from x = to x = 7, we find 

1 n 

-q od = (E od - E b ) / Xod{x)Xb{x)dx. 



From (A.27)-(A.29), we see that 

p b (3 = 7T - 6 b = 71 and 6 b = 6^ = 

Likewise, we can also show that 

PodP = n - 9 od = n and 

Thus, q od = q b = W, and the integral in (A. 99) is 

n 

Xod(x)xb(x)dx = / Xod(x)xb(x)dx 



IT 



Wf3 



od — "co — 



7T 

wp' 



sinh q od a cosh q b a p 
sinp od/ 3 sinp 6 /? 

-(3 = 



I sin 2 poo ( —x + 7) drr 



40 O A 2 " 8tt 2 
Since g od = W^, we derive from (A. 91) 

_ 1 2vr 2 



(A98) 



(A99) 



(A100) 



(A101) 



(A102) 



(A103) 



in agreement with (A. 86). 

We are now ready to introduce the two-level model. We shall approximate the Hamiltonian 

T + V(x), T + V a (x) and T + V b (x) of (A.7)-(A.9) by the following three 2x2 matrices: 



h = E 00 + 



-A 



- A 




(A104) 



2 

and 



h a = E 00 + ^ 2 + ( _° A A ) (A105) 



h = E oo +( x A j, (A106) 
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with ip, Xa and Xb as their respective groundstates which satisfy 

hlp = Elp, h a Xa = E a Xa 

and h bXb = E bXb - (A107) 

The negative sign in the off-diagonal matrix element —A in (A.104)-(A.106) is chosen to make 



Xa Xb 



V2 { 1 



(A108) 



simulating the evenness of Xa{x) and Xb(x). Likewise, the analog of Xod is the excited state 
of h b , with 



Xod V2 { "I 



and 



hbXod — E od Xod- 



It is straightforward to verify that 



where 



sin 2£ 



4A 



sin£ 
cos£ 



and cos 2^ 



V( 4A ) 2 + ^ 



(A109) 
(A110) 

(Alll) 
(A112) 



and 



E, 



When A << ^ , we have 



in agreement with (A. 88). 



E = Eoo + — 
E a = E^ + ^ — A 

Eh = Eno — A 



A2 1 T 



od - c 'oo 



^ = -&oo T" 

n 2 



(A113) 



(A114) 
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Next, we wish to examine the relation between the two-level model and the soluble square- 

2 

well example when A is O(^). Assume, instead of (A. 76), 

{Way 1 « i, (wpy 1 « 1, e- 2q ~ a « 1 

but (A115) 

e -2 qoo 5 _ Q^y 

Hence, in the square-well example, (A. 83), 

Eoo -E = ^Le- 2 ^ a+8 \ 



W(3 3 

and (A.86), 



Wf3 3 

remain valid; on the other hand, (A. 54) and (A. 78) now lead to 



— T T 7" 03 ' 



L 2 ^ \( e 2q °° s - e- 2qocS ) = 2Asmh2 ?0O <5. (A116) 

Thus, the above expressions for i^oo — _E and A give 

E^-E = A e~ 2qac5 = A (cosh 2^^^ - sinh 2 qoo 5) 
= \(\J 1 + sinh 2 2goo5 — sinh 2q O0 5). 

Together with (A. 116), this shows that the soluble square-well example yields 



E^ — E 



in agreement with (A. 113) given by the two-level model. 

In both the square-well problem and the simple two-level model, we can also examine the 

2 

limit, when A >> In that case, (A. 113) gives 



E = E 00 -X + ^--^ + 



4 32 A 



which leads to 

E = Eoo - A = E b when /j 2 = 0, 

in agreement with the exact square-well solution. Furthermore, if we include the first order 
correction in 0(/i 2 ), (A. 115) gives 

E = E OQ + ^-X = E b + X -{E a - E b ) + O(^). (A117) 
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As we shall discuss, for the exact square-well solution, (A. 117) is also valid. Thus, the simple 
two-level formula (A. 113) may serve as an approximate formula for the exact square-well so- 
lution over the entire range of ^/i 2 /A. 

A. 3 Square-well Example (Cont.) 

We return to the soluble square-well example discussed in Section A.l. As before, ijj(x) 
is the groundstate of T + V(x) with energy E, which is determined by the Schroedinger 
equation (A. 7). Likewise, x( x ) ,s the trial function given by (A. 57); i.e., the groundstate of 
T + V(x) + w(x) with eigenvalue E = E a , in accordance with (A.58)-(A.60). From (A. 59) 
and (A. 65), we see that the energy difference 



satisfies 



where 



and 



S = E a -E (A118) 
M — P x(x)ip(x)dx (A120) 

JO 
,0 

N= / x (x)4>(x)dx. (A121) 



Before we discuss the iterative sequence {£„} that approaches S, as n — > oo, it may be 
instructive to verify (A. 119) by evaluating the integrals (A. 120) and (A. 121) directly. Choose 
the normalization convention of ip and \ so that at x = 7 



X'il) 

From (A.10)-(A.12) and (A. 57) we write 



(A122) 



sin k (—x + 7) a < x < 7 

^ x) = ^f^a\ J (^H^)cosh q (x-5) for -a<x<a (A123) 



k \ smk a (3 




sin kB cosh q(a+S) \ / . \ . . 

sin P gcosMa-a) J sm ^ + 7) -7 < x < -a 



a < x < 7 



cosh q a x 
cosh q b x 



for 



< x < a 
— a < x < 



(A124) 



— 7 < rr < —a. 
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By directly evaluating the integral / x{x)ip{x)dx , we can readily verify that for 7 > x > 

\{k 2 - hi) £ X (vMv)dv = \{x{xW{x) - 4>(x) X '(x)) (A125) 

and for — 7 < x < 0, 

1 f x 1 

2 (P 2 - Pb) J ^ X(y)i'(y)dy = --(x(x)ip'(x) - ip(x)x'(x)). (A126) 

Both relations can also be inferred from the Schroedinger equations (A. 7) and (A. 58). Setting 
x = and taking the sum (A.125)+(A.126), we derive 

1 -(k 2 -kl)M+^- P l)N = 

which, on account of (A.13)-(A.18), lead to the expression for the energy shift 8, in agreement 
with (A.119). 

Next, we proceed to verify directly that f(x) = vp(x)/x(x) satisfies the integral equation 
(A. 68). With the normalization choice (A. 122), we find at x — 7, since ^(7) = x{l) — 0> 

m = ^ • (A.127) 

which gives the constant in the integral equation. The same equation (A. 68) can also be cast 
in an equivalent form: 



f(x) = 1 + f x-\x){x\G\z)x{z){w{z) - S)f(z)dz (A128) 
j -7 



H 
-7 

where (x\G\z) is the Green's function that satisfies 



(T + V(x) + w(x) - E )(x\G\z) = S(x - z) 
and (A129) 

(x\G\z) = for i>z. 

For x < z, (x\G\z) is given by 

(x\G\z) = -2( X (x)x(z) - x(x)x(z)) (A130) 

where 

= X (x) / (A131) 
is the irregular solution of the same Schroedinger equation (A. 58), satisfied by x( x )- I.e., 

(T + V(x) + w(x))x(x) = E x(x). (A132) 
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Consequently, over the entire range — 7 < x < 7 

X'(x)x(x) - x'(x)x(x) = 1. 
According to (A. 11), (A. 12) and (A. 57), we have 

' A sin ka (_ x + 7 ) + _L (jg^L) cos ^(-^ + 7 ) 

- sinh o a a; 

Y sinh g^a; for 

. - B <**Pb(x + 7) - i (Jf^) cosp b (x + 7) 
where A and B are constants given by 

— J_ ( sinh q a a \ 1_ ( cos k a (3 \ 

q a V sinfc a /3 J k a \cosh q a a) 



(A133) 



a < rr < 7 
< x < a 

— a < x < 

— 7 < rr < —a 



(A134) 



and 



_1_ / sinh <?t,a \ 1_ / cos p b f3 \ 

~ q b \ sin p b /3 / p 6 Vcoshg^aJ 



(A135) 



Since in (A. 128), there are only single integrations of the products x( z )^>( z ) and xl- 2 )^!- 2 ). 
one can readily verify that f(x) satisfies the integral equation, and therefore also its equivalent 
form (A.68). 

A. 4 The Iterative Sequence 

The integral equation (A.68), or its equivalent form (A. 128), will now be solved iteratively 
by introducing 

ipn(x) = x(x)fn(x). (A 136) 

As in (4.87)-(4.89), f n (x) and its associated energy £ n are determined by 



and 



When n = 0, we set 
Introduce 

and 



f n (x) = 1 -2 !\- 2 {y)dy f x 2 {z){w{z)-£ n )f n ^{z)dz 



I" X 2 (z)(w(z)-S n )f n - 1 (z)dz = 0. 

J-y 

fo(x) = 1. 

Mn = ^ X 2 {x)f n {x)dx 
JO 

Nn= X 2 (x)f n (x)dx. 
J-7 



(A137) 

(A138) 

(A139) 
(A140) 

(A141) 
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From (A. 59) and 
we derive 



E — E — S 



N _i 



and 

_ M n _ x E a + Nn-iEb 

For n = 1, we have from (A.139)-(A.141), 

if 1 . , „ cosh 2 g a' 1 . ] 

M ° = 2 r + ^ smh2 «° a + ( " - 2K sm 2 ^> } • 

if 1 . . n cosh 2 g b o; 1 . 1 

and 

i?i = H — Eh)- 

1 6 M + No V a b! 

For small fj 2 , since i? a — E b and M — N are both 0(fi 2 ), we find 

E 1 ^E b + ^(E a -E b ) + 0(f J , 4 ) 

in agreement with (A. 117), given by the simple two-level formula. 

Next, we examine the integration for f n (x). Consider first the region 

a < x < 7; 

(A. 137) can be written as 

f n (x) = l + 2S n \\-\y)dy \\\z)U-x{z)dz. 

Jx Jy 

Introduce 

i = k a (-x + -f) 

and 

When n = 0, we set 

MO = sin £- 
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From (A. 149), or more conveniently by using given by (A. 130), one can readily verify 

that, for a < x < 7, 



vi(0 



1 1 
(1+2 e sin ^ — 2 61 ^ C0S ^' 



f 1 3 1 I 13 

|(1 + ^ e 2 + g eie 2 ) - - eie 2 £ 2 jsin£- (- e 2 + - eie 2 )£cos£, 

|(1 + ^ £3 + g £263 + eie 2 e 3 ) - ^(e 2 e 3 + eie 2 e 3 ) £ 2 J sin£ 

f A 3 5 x , 1 x$l 

+ e 3 + g e 2 e 3 + 7^ eie 2 e 3 ) £ + — eie 2 e 3 £ \ cos£, 



(A154) 



etc. These solutions can also be readily derived by directly using the differential equation 
satisfied by ip n (x) = x(x)f n (x): 



(T + V(x) + w(x) - E a )ij n {x) = (w(x) - 4)^n-i(^), 
where in accordance with (A. 14), E a = |(/i 2 + k 2 a ). For a < x < 7, we have 



(A155) 



V(x) 



w(x) = 



and therefore 



Introduce S n (£) and C n (£) to be polynomials in £, with 

v n (0 = (f[e m ) {S n (0 sine + C n (0 cos£} . 
1 

From (A.152) and (A.156)-(A.157), we find 

5 n (0-2C n (0 = S n _i(0 



(A156) 



(A157) 



and 



(A158) 



C„(0 + 25„(0 = C n _x(0, 



where the dot denotes so that 5 n = etc. At x — 7, we have £ = 0, / n (7) = 
^n(7)/x'(7) = 1 and therefore 



S n (0) + C n (0)= He' 1 . 



(A159) 



m=l 



For n = 0, So(£) = 1 and C (£) = 0. Therefore, for n—1, (A. 158) becomes 



Si - 2Ci 
Ci + 251 



1 




(A160) 
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Assuming Si and C\ to be both polynomials of £, we can readily verify that Si is a constant 
and C\ is proportional to £. Using (A. 160) and the boundary condition (A. 159), we can 
establish the first equation in (A. 154), and likewise the other equations for n > 1. 

To understand the structure of t>i(£). ^(O- ^(O; ' ' we mav turn t0 tne e xac t solution 
ip(x) given by (A. 123). In analogy to (A. 152), we define through 

«.) S ^.<0. (^.161) 

Thus, for a < x < 7, 

u(f ) = ^ sin + 7). (A162) 
From (A.13)-(A.14) and (A.118), we have 

t = E a -E = \ (k 2 a -k 2 ). (A163) 

In terms of 

2f k 2 

^|=l-|. (AW) 

we write 

v(0 = -j=L= sm^VT^e) (A165) 
VI - e 

with £ given by (A. 150), as before. It is straightforward to expand as a power series in e: 

+ {-(^ + ^ 2 + ^ + le 3 e 3 } cosC + 0(e 4 ) (A166) 

To compare the above series with i> n (£) of (A. 154), we can neglect 0(e n+1 ) in (A. 166). The 
replacements of all linear e-terms by e„, e 2 -terms by e n _ie„, e 3 -terms by e n _ 2 e„_ie n , etc. lead 
from (A. 166) to i> n (£)- ' s °f interest to note that the expansion (A. 166) of v (£) in power 
of e has a radius of convergence 

|e| < 1. (A167) 

On the other hand, the iterative sequence ' s always convergent, on account of the 

Hierarchy Theorem. The main difference between (A. 154) and (A. 166) is that in (A. 154) each 
iterative e„ is determined by the fraction (A. 142). 
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region 


ei(x) 




ipi(x) 


I 

II 
III 

IV 


|(£ a -#i) 
^(E b - Ei) 

§z( E b - E x ) 


k a (-x + 'j) 
q a x 
-q b x 

Pb(x + 7) 


^{(l + ^sme-l^cose} 
(kh + |ei£) sinh£ + p n cosh£ 
(/tin + |eiO sinhf + p m cosh£ 
^ {Kv + i Cl ) sine -fa cose} 



Table 2. The n = 1 iterative solution in the four regions: 

I (a < x < 7), II (0 < x < a), III (-a < x < 0) and IV (-7 < x < -a). 
The constants E\ = E a — Si, k\\, Km, kiv, pu and pm are given by (A.175)-(A.177). 

In a similar way, we can derive ip n (x) in other regions, —a < x < a and — 7 < rr < —a. 
The results for n = 1 are given in Table 2. The functions ei(a;) and £(x) are discontinuous 
from region to region. The constants K n and p n are determined by requiring ipi(x) and ^i(x) 
to be continuous at x = a. In region I, when x = a+, we have 

cosh q a a 



and 



^i(oH-) 



^i(a+) 



sin fc a /3 



(1 + ^ei) sin k a (3 - ^eik a /3 cos fc a /?J 



-/c a CQShga " jcos fc a /3 + ^eifc a/ 9 sin fc a /3 

Sill rZaP ^ ^ 



where the constant 
with 

In region II, when x = a— 
and 



(E a -E 1 ) = e 1 (x) in I 



Ei = E a — Ei. 



(«n + 2 e ii?a«) sinh g a a + p n cosh g a a 



^i(a-) = g a |(/c n + -e u q a a) coshg a a + (p n + 7^11) sinhg a a| 



where the constant 



en = —(E a - Ei) =6i(x) in II. 
9S 



The constants k\\ and p n are determined by 

ijji(a—) = ipi(a+) and 
Likewise, the constants Km and pm are given by 

^i(O-) =^i(0+) and 
and the constants k iv and E x are determined by 

ijji(-a-) — ijji(-a+) and ip[(-a-) = ip[(-a+). 



^i(O-) = vi(o+), 
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